Hierarchical Reference Theory of critical fluids in disordered porous media 
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We consider the equilibrium behavior of fluids imbibed in disordered mesoporous media, including 
their gas-liquid critical point when present. Our starting points are on the one hand a description of 
the fluid/solid-matrix system as a quenched- annealed mixture and on the other hand the Hierarchical 
Reference Theory (HRT) developed by A. Parola and L. Reatto to cope with density fluctuations 
on all length scales. The formalism combines liquid-state statistical mechanics and the theory of 
systems in the presence of quenched disorder. A straightforward implementation of the HRT to 
the quenched-annealed mixture is shown to lead to unsatisfactory results, while indicating that the 
critical behavior of the system is in the same universality class as that of the random-field Ising 
model. After a detour via the field-theoretical renormalization group approach of the latter model, 
we finally lay out the foundations for a proper HRT of fluids in a disordered porous material. 

PACS numbers: 11. 10. Hi, 75.40.Cx 



I. INTRODUCTION 

The HRT developed by A. Parola and L. Reattoi^.2. 
combines liquid-state statistical mechanics and nonper- 
turbative renormalization group (RG) formalism. As 
such, it has proven a powerful tool for describing the 
equilibrium structure and the phase behavior of simple 
fluids and fluid mixtures in a way that properly incor- 
porates the long wavelength fluctuations present in the 
vicinity of the critical point (s)i"—. 

The purpose of this article is to show how the HRT 
can be generalized to deal with fluids in disordered meso- 
porous media. The critical behavior of such systems is 
indeed an unsettled question. It was suggested by de 
Gennes and Brochard- that fluids imbibed by very di- 
lute disordered porous media like silica gels represent an 
experimental realization of the random-field Ising model 
(RFIM). The rationale is as follows: as the solid attracts 
preferentially one phase of the fluid (usually the liquid 
for a single-component fluid and one of the species for a 
mixture) , there is a shift of density or concentration away 
from the solid; this can be described as the result of a 
perturbation to the chemical potential which, due to the 
disordered structure of the gel, is random in space. At a 
coarse-grained level, the fluid near criticality in the pres- 
ence of a gel can then be modeled by a lattice gas with a 
random chemical potential which, provided the correla- 
tion length of the fluid is much larger than that charac- 
terizing the structure of the gel, can be taken as uncor- 
related from site to site. Going to the Ising-spin repre- 
sentation of the lattice gas thus leads to a ferromagnetic 
Ising model in a random magnetic field. The random- 
field Ising model (RFIM)^ is an archetype of statistical- 
mechanical model in the presence of "quenched disorder" . 
The quenched disorder in the fluid case is associated with 
the structure of the solid matrix of the porous medium. 

It has however proven extremely difficult to test the de 
Gennes-Brochard suggestion in experiments and in com- 
puter simulations. This is due to strong crossover effects 



to bulk-like behavior when the effect of the quenched dis- 
order is very weak^ and to extremely slow dynamics in 
generali^— . Actually the vast majority of experimen- 
tal studies appear to be out of equilibrium and charac- 
terized by the presence of hysteresis effects. Such an 
hysteresis is routinely found in adsorption experiments^^ 
where the filling of the porous medium by the fluid 
(the "capillary condensation") and the draining of the 
medium (the "capillary evaporation") do not take place 
at the same chemical potential. Phase transition could 
still be observed in very open porous materials such as 
light aerogelsiS., but these are out-of-equilibrium phase 
transitions!^'!^. 

The lack of strong experimental evidence showing the 
equivalence between the critical behavior of fluids in dis- 
ordered mesoporous media and that of the RFIM could 
come from several features that are missing in the sim- 
ple RFIM: the fact that the porous material imposes a 
correlated disorder on the fluid or that the distribution 
of random chemical potentials is strongly asymmetrioiS, 
the effect of confinement and/or of dilution, etc. On the 
other hand, a recent series of Monte Carlo studies of sev- 
eral models of fiuids in disordered medians has concluded 
in favor of the de Gennes-Brochard conjecture. 

The HRT appears as a promising candidate framework 
to tackle the behavior of fluids in disordered mesoporous 
materials. It is indeed able to provide a nonperturbative 
implementation of the RG in the appropriate regime near 
the critical point. It also allows one to start with, and 
maintain, a microscopic description of the system, in- 
cluding in the present case confinement and/or dilution, 
correlated asymmetric disorder, and wetting phenomena. 
One should finally add that the critical behavior of the 
RFIM has itself been under debate since its introduction 
in the mid-seventies^ and that the standard perturba- 
tive RG approaches have been shown to be inadequate to 
describe the critical behavior of the system^. Nonpertur- 
bative RG methods then offer a powerful alternative^"—. 

The outline of this paper is as follows. In section [H] 



we recall the basics of the HRT formalism by using a 
functional formalism and notations that are useful to 
make the connection with the field-theoretical RG of dis- 
ordered systems. Section [IIII is a summary on the micro- 
scopic description of fluids in disordered porous media 
known as the "quenched-annealed mixture" and of the 
so-called "replica formalism" borrowed from spin-glass 
theory. The next section, section IIVI is a first, which 
after the facts can be dubbed as "naive", attempt to 
generalize the HRT to the quenched-annealed mixture. 
We derive the HRT hierarchy and its asymptotic form 
in the vicinity of the critical point and at long wave- 
lengths. The very same procedure used for the RFIM 
leads to flow equations that are formally similar to those 
of the quenched-annealed mixture, giving credit to the 
equivalence of the critical behaviors of the two systems. 
We also present a flrst application in which we analyze 
the asymptotic RG flow equations for the fluid within an 
Ornstein-Zernike approximation. The results display in- 
consistent and unsatisfactory features. Section IVl then is 
a discussion of the flaws of the naive HRT and a short 
review of the lessons that can be drawn from a recent non- 
perturbative RG study of the RFIM. In section lVTl we lay 
out the foundations of a proper HRT treatment of fluids 
in disordered mesoporous materials, deriving the hierar- 
chy of HRT flow equations, their asymptotic form, and 
proposing approximation schemes. We flnally conclude 
in section IVIII and we provide some additional technical 
details in two appendices. 



II. HRT IN BRIEF 

The HRT incorporates in the realm of liquid-state the- 
ory a proper account of the long wavelength fluctuations 
that play a central role in criticality This is done by pro- 
gressively including longer wavelength (shorter wavevec- 
tor) components of the attractive part of the interac- 
tions in the computation of the equilibrium properties 
of a fluid. In the original HRT formulationii^, this was 
achieved by means of a "sharp infrared cutoff" that pre- 
vents integration over all wavevectors strictly less than 
a running value k. Since then, it has been reformulated 
in a more general framework that can accommodate any 
smooth form of infrared (IR) cutoff procedure^i^. To 
make an easier contact with the fleld-theoretical RG, 
most notably that of disordered systems, we use the more 
compact functional formalism and change notations from 
the conventional HRT ones; the correspondence with the 
latter is given in the text and recalled in table HI 

Consider an atomic liquid at equilibrium in the pres- 
ence of an inhomogeneous chemical potential fi{x), where 
X represent the coordinates of a point in d-dimensional 
space. The atoms interact via a pair potential that can 
be decomposed as v{\x — y\) — vr{\x — y\) + w{\x — y\), 
where v^ is a steep repulsive interaction (often modeled 
as a hard core) and w is a longer-ranged attractive in- 
teraction that can trigger a liquid-gas transition with a 



critical point. The grand partition function of the fluid 
reads 



S[/i]=X!Tn/ ••■ / exp{- I3[VbXxi,...,xn) + 

jY -'* • Jxi Jxn 



1 



w{\x ~ y\)p{x)p{y) - / ii{x)p{x)\] 



(1) 



where Vr is the sum of the pair interactions w^ among all 
the fluid particles [Vr = (1/2) /^ /^ vr{\x - y\)p{x)p{y)], 

P(^) = ^i=i^'''^K^ ~ ^i) is the microscopic density of 
particles at point x, J = J d'^x, and p{x) may be cor- 
rected if necessary by a self-energy term — w(0)/2-'^. 

The idea behind the exact or nonperturbative RGi^"— 
is to progressively include the effect of fluctuations (more 
specifically here, the fiuctuations generated by the attrac- 
tive part of the interactions) by integrating over longer 
and longer wavelengths or, alternatively, smaller and 
smaller wavevectors. This is achieved by introducing an 
IR regulator that suppresses the integration over fluctu- 
ations of wavevector (momentum) smaller than a cutoff 
k. At the beginning of the RG flow, fc = A represents 
a microscopic ("ultra-violet") scale whose inverse in the 
present case is of the order of the typical interatomic dis- 
tance: at this scale, virtually no fluctuations are taken 
into account, as in a mean-fleld approximation. On the 
other hand, when fc = 0, all fluctuations have been in- 
cluded and one obtains the exact theory. 

The HRT combines this RG feature with the require- 
ment of a proper description of the (noncritical) physics 
mainly induced by the short-ranged repulsive forces in 
liquids; the IR regulator only modifies the attractive in- 
teraction between particles. Let us then define in Fourier 
space 



Wkiq) = w{q)Rk{q) 



(2) 



where Rk{q) is an IR regulator that goes to zero when 
k ^ Q and is equal to 1 when fc = A; by construction 
then, Wk=A{q) = w{q) and Wk=t){q) = 0. The regulator 
is also chosen so that it goes to zero as g S> fc and to a 
fc-dependent constant (a "mass" in field-theoretical lan- 
guage) when q <^ k, which will be helpful to limit the 
contribution of the low-g modes in the grand partition 
function without affecting that of the high-g modes. 

Note the changes of notation with respect to the orig- 
inal HRT literature: the IR cutoff fc and the momen- 
tum q are denoted by Q and fc, respectively, in the lat- 
ter and the Wk defined here corresponds to w — wq in 
Parola-Reatto's work^"^. In addition, we have followed 
Caillol^^ in the way to introduce the regulator. (Parola 
and ReattoS, have a slightly different procedure that uses 
the functional form of w(r) to define the cutoff.) The rea- 
sons for our choice of notations comes from their closer 
connection to what was done in the field-theoretical non- 
perturbative RG, more specifically when applied to the 
theory of disordered systems. 

One defines at the scale fc the grand partition function 
of an equilibrium fiuid whose atoms interact through the 



pair interaction vp 



Wk as 






{xi,...,xn) + 



[w{\x - y\) - Wk{\x - y\)\p{x)p{y) 
H{x)p{x)\], 



(3) 

so that, as a consequence of the properties of the IR reg- 
ulator, Sfc is equal to the grand partition function of the 
reference system, Sfc=A[/^] — Si^[/i], at the microscopic 
scale and is equal to the exact grand partition function 
of the fully interacting system, Sfc^oiA*] = ■^'M with S[/i] 
given in Eq. ([T]), at the end of the flow. 

From the grand partition function at the scale fc, one 
defines the "effective average action" in field-theoretical 
language^ (which, up to a factor /3, is the Helmholtz free- 
energy functional in liquid-state theoretical language) 
through a Legendre transform followed by the addition of 
a mean-field-like contribution that corrects for the pres- 
ence of the regulator: 



r.M 



+ 



WM+M t^{x)p{x) 

J X 

'Wk{\x-y\)p{x)p{y), 



(4) 



where Wfc[/i] = logE!fe[/x] is the generating functional 
of the connected many-body correlation (or Green's) 
functions (up to a factor — /3, Wk is a grand-potential 
functional in liquid-state theoretical language), p{x) =< 
p{x) > is the average (inhomogeneous) density field ob- 
tained from the Legendre transform as 



p{x) 



5WM 
I36p{x) ■ 



(5) 



(If necessary, one can take account of the self-interaction 
terms by replacing p{x) by p,{x) — Wk{0)/2.) The effec- 
tive average action (modified Helmholtz free-energy func- 
tional) Ffc is the generating functional of the "proper" 
or "1-particle irreducible" (IPI) vertices^l, which are es- 
sentially the direct correlation functions of liquid-state 
theory. At the microscopic scale fc = A, it reduces to a 
mean-field description of the attractive interactions. 



Tk=A[p] =Tb.[p] 



w{\x-y\)p{x)p{y) (6) 



where Tj^ [p] is the effective action of the reference system. 
[Note that one could choose A — ?> oo, but the fiuctuations 
generated by the attractive interactions with wavevectors 
much larger than 2TT/a, where a is the typical interatomic 
distance, essentially do not renormalize eq. (JH).] At the 
macroscopic scale fc = 0, F^ is equal to the effective action 
of the fully interacting system. 



Ffc=o[p] = -\og=.[p] +(3 p{x)p{x) 



(7) 



It is worth stressing that due to the addition of the last 
term in Eq. (HI, F^ is not exactly the Legendre transform 
of Wk and therefore need not be convex except when 
fc->0. 

The 2-point (connected) correlation function is ob- 
tained from VFfe[/i] as 



Wlt\x,y) 



5^Wk[p] 
(3Sp{x)pSp{y) ■ 



(8) 



It is related to the standard 2-point total correlation func- 
tion hk{x,y) of the inhomogeneous hquid at scale k, in 
which the atoms interact with the modified pair potential 
vr + w — Wk, according to 

Wlf\x,y) = pix)S^''Hx - y) + p{x)p{y)hk{x,y). (9) 

(M^(2) is often denoted as F or F^^) in the HRT 
literaturei"— .) For a uniform density field p{x) = p, its 
Fourier transform, Wl^\q,q') = {2TTY5^'^\q+q')wl^\q), 
is then connected to the static structure factor Sk{q) 
through W^^^q) = pSk{q)- Similarly, the 2-point "IPI 
vertex" 



rf(x,y) = 



5^Vk[p] 
5p{x)Sp{y) 



(10) 



is related to the standard Ornstein-Zernike (OZ) direct 
correlation function Ck{x,y) of the liquid characterized 
by the modified interaction at scale k according to 

Ff(x,y)-/3^fc(|x-y|) = ^^^^^p^-c,(x,y). (11) 

The 2-point functions VF^ {x, y) and F]^ '{x, y) are hnked 
through the (modified) Legendre transform as 



.(2) 
k 



Wlf\x,y) ^[T['^ ~ pwk]-\x,y), 



fl2) 



where the inversion is understood in terms of operators. 
The above equation when specified to a uniform density 
(and equivalently uniform chemical potential) is nothing 
but the familiar OZ equation for the (modified) liquid at 
scale fc: 1 -I- phk{q) = [1 - pck{q)]~^ ■ 

Upon decreasing the infrared cutoff fc, the effective av- 
erage action (modified Helmholtz free-energy functional 
at the scale fc) F^. evolves according to the following exact 
"flow" equation: 



dtTk[p\ = \ 



dt(i>ki\x - y\) Vf\p\+<t>k {x,y), 



X Jy 



where t — log(fc/A) and we have (conventionally) intro- 
duced the notation (j)k — —/3wk, which is a positive quan- 
tity as Wk is an attractive interaction. By differentiating 
Eq. ([T^ with respect to the density field and evaluat- 
ing the resulting expressions for a uniform density field 
p{x) = p, one generates an exact hierarchy of coupled 
flow equations for the IPI vertices (modified direct cor- 
relation functions), which is the more familiar form of the 



HRT. For instance, the first equation for the Helmholtz 
free-energy density Ak{p) = rk{p)/V (or "eflFective aver- 
age potential" in field-theoretical language) reads 



dtAkip) 



pdt(j)k{q) 



1 



2./g l-pcfc(g;p)' 



(14) 



where J = <P/{2ti)''- and Ck{q; p) is the direct OZ corre- 
lation function that is related to the 2-point IPI vertex 

hyl/p-Ck{q;p)=ri^\q;p) + Mq) [see Eq. (HI])]. Note 
that we have defined Ak as the opposite of the defini- 
tion chosen by Parola and Reattoi"— : in our notation Ak 
is directly /3 times the modified Helmholtz free-energy 
density at scale k. 

The above functional integro-differential equation for 
the RG flow of the effective average action, Eq. p^ . is 
exact but of course impossible to solve in general. Ap- 
proximation schemes must be devised. In field theory, 
nonperturbative approximations are formulated directly 
at the level of the generating functional^i^S: an ansatz 
is chosen for the effective average action and consistency 
is automatically enforced for the relations between the 
derivatives of the latter and the IPI vertices. In the 
HRT of fluids, it is not as easy to work at the level of 
the generating functional, here the Helmholtz free-energy 
functional; approximations are proposed as truncations 
of the hierarchy of coupled equations for the IPI vertices 
for a uniform density field. The common implementation 
corresponds to formulating an ansatz for the 2-point IPI 
vertex function, here the direct pair correlation function. 
One has then to be careful to properly enforce the sum 
rules that relate the different levels of the hierarchy. For 
instance, the 2-point direct correlation function must sat- 
isfy the "compressibility sum rule" that relates it to the 
Helmholtz free energy, namely, 

i - c, (g = 0; p) - 0, (g = 0) = rf {q = 0; p) = ^^0^ , 

where we have used the definitions and notations intro- 
duced above. 

A typical closure of the HRT hierarchy can be written 
asi-^ 

- - ck{q: p) - Mq) = ^^44^ + ^^'^^'i' ^)' (16) 

p op^ 

where an approximation is formulated for 6T.k{q', p) that 
satisfies 5Yjk[q — 0; p) ^ 0. Actually, the theory being 
regularized for wavevectors q ^ k, it is expected that a 
relatively simple ansatz for Ck{q',p), or equivalently for 
6Y.k{q',p), already captures most of the physics. All of 
the proposed approximations have been of "OZ type" as 
they assume that the ^-dependence of 5T.k{q', p) is ana- 
lytic at small g's and therefore starts as q^. The simplest 
nontrivial such approximations are of RPA form with e.g. 



6T,k{q; p)=- [cR{q\ p) - CR{q = 0; p)] 
- [0(g) - cj^iq = 0)] , 



6Y.k{q]p) = -[cfl,(g;p) 
>(g)~0(g^O) 
0(9 = 0) 



CR{q^ 0;p)] 

d'Akip) , , n ^ 1' 

— ^-^ h CR{q = 0; p) - - 

op^ p 



(18) 

The latter expression comes from the OZ-like closure 
Ckiq; p) + 4>kiq) = CR{q; p) -|- Xk{p)(l>{q), where the T and 
p dependent factor Afe is adjusted so that Eq. (fTSJ) is sat- 
isfied. Improved approximations based on ORPA, MSA, 
etc, which allow one to enforce the "core" condition when 
the reference pair potential is a hard-sphere one, have 
also been implemented^"— . 

The main drawback of all these OZ-like approxima- 
tions is that they do not properly describe the singular 
wave- vector dependence of the correlation function at the 
critical point, i.e. the fact that when fc = 0, dT,k{q) starts 
as q'^~'^ instead of q^, where 77 > is a critical exponent 
called the "anomalous dimension" of the field. [(5Sfc(g) 
is analytic for q ^ k but the singularity emerges when 
g ~ fc — >■ 0.] This weakness is not too significant when 
describing the critical behavior of pure fluids or mixtures, 
as rj ~ 0.03 is then very small, but it is much more serious 
when considering the critical behavior of a fluid in a dis- 
ordered porous medium. This will be further discussed 
below. 

Finally, we close this brief review on the HRT by con- 
sidering the sharp cutoff originally proposed. In this 
case, the regulator is then a Heaviside step function, 
Rk{q) oc H{k—q), so that Wk = when q > k and Wk = w 
when q < k. Its derivative dtRk{q), which appears in the 
exact flows equation Eq. (1131) . is then proportional to a 
delta function, S{k — q). Some care is needed when han- 
dling the "loop integral" appearing in the right-hand side 
of the flow equations as it involves a combination of delta 
and Heaviside functions [see for instance Eq. ([Ti| ]. This 
point is detailed in Rcfs. [32, [25. After a proper account 
of this problem, Eq. (J14p can be rewritten as 



dtAkip) ^^k'^vdlog 



0(fc) 



rf(A:;p), 



(19) 



where v^'^ = 2'^+^Tr'^/'^T{d/2) and (p{k) = ~l3w{q = fc). 
The above flow equation is identical to that derived by 
Parola and Reatto (with the change of notations Q ^)- k, 



A, 



Q 



^Ak, Cq^T 



(2) 



and • 



' ~" '/'/c)- The same is 



(17) 



true for the higher orders of the HRT hierarchy. 



III. QUENCHED-ANNEALED MIXTURES, 

REPLICA FORMALISM, AND LIQUID-STATE 

THEORY 

The first step in the description of the behavior of fiuids 
adsorbed in disordered mesoporous materials is to pro- 
vide a realistic model of the latter. The approach that 
has been used in the past two decades, both in computer 
simulations and in statistical-mechanical theories^i^, is 



field — theoretical RG 


HRT 


IR cutoff k 


IR cutoff Q 


Momentum q 


wavevector k 


Effective average action Fk 


luodifled Helmholtz free-energy functional 


Effective average potential rk{4>)/V — Ak 


modified Helmholtz free-energy density —Aq{p) 


Proper (IPI) vertices: at the 2-point level, f). 


Modified direct correlation functions: at the pair level, Cq 


Green's (correlation) functions: at the 2-point level, Wl 


(Connected) correlation functions: at the pair level, Fq 


Regulator for fluctuations: (j>k{q) = —l3Rk{q)w{q) 


Modified interaction: (j){k) — 4>Q(k) = "I3[w{k) — WQ{k)\ 


Legendre transform relation: WJ. = {Tj, + (pk)^^ 


Ornstein-Zernike equation: Fq — {—Cq + 4> ^ 4>q)^^ 



TABLE I: Table of equivalence between the field-theoretical nonperturbative RG and the HRT for fiuids. 



to treat the solid phase of the porous material (the "ma- 
trix" ) as a collection of particles in some predefined mi- 
crostructure which is assumed to be rigid'^^''^^. There- 
fore, the solid is simply vievifed as a configuration of 
quenched particles, sampled from a given probability dis- 
tribution. This distribution can be obtained either by 
modelling the process of formation of the porous mate- 
rial and considering the resulting configurations (as done 
for instance for the porous glass Vycor formed by spin- 
odal decomposition^ or for the silica aerogel modelled 
by a cluster-cluster aggregation'^^''^^) or by taking some 
chosen equilibrium distribution of particles (hard spheres 
or ideal gas). 

The fluid establishes itself in a state of thermal equi- 
librium in the presence of, and in interaction vifith, 
this quenched structure. Such a "quenched-annealed 
mixture"—'^ is macroscopically homogeneous, but the 
external field exerted by the matrix on the fluid varies 
on micro- and meso-scopic scales, thus making the fluid 
strongly inhomogeneous on such scales. This ham- 
pers the direct application of many common methods 
of liquid-state theory. For instance, density-functional 
schemes, which have proven very useful for fluids in 
the presence of a substrate of simple geometry^, are 
intractable in general for a realistic disordered porous 
medium (the required numerical resolution becomes only 
barely possible for simplified lattice-gas model o^^'^^i"^° ). 

It was shown^^i^ that the difficulty can be, at least 
formally, circumvented by applying the "replica method" 
developed in the theory of spin glasses^. The study of 
the original quenched-annealed mixture is then replaced 
by that of a (fictitious) homogeneous mixture of n -I- 1 
components at equilibrium: one component is formed by 
the matrix particles now considered as fully annealed, 
and the others are n identical replicas of the fluid. The 
matrix interacts with all fluid replicas, but the latter do 
not interact with each other. For this (n -I- l)-component 
mixture, one can use the whole machinery of liquid-state 
statistical mechanics and, at the end of the manipula- 
tions and calculations, take the n — ?■ limit with an 
appropriate analytic continuation. 

Consider for instance the grand potential Vtf (/i) of the 
fluid in equilibrium at a chemical potential /i within a 



porous material [in the fleld-theoretical language used in 
the above section, —f3ilf{fi) — W/(/i)]. It is obtained as 



-PQfip) 



Po{{yt}No)^og:if{n;{yi}No) 



Vi^o 



logS/(Ai;{2/i}Aro), 



(20) 

where {yi}No denotes a configuration of iVo matrix parti- 
cles, 'PodyiJA'o) is ^^6 probability of finding this config- 
uration, and the overline in the second line is a compact 
way to express the average over the "quenched disorder" , 
i.e., the frozen-matrix configurations; 'E,f(fi;{yi}j^g) is 
the matrix-dependent grand partition function which is 
defined as 



S/(M;{2/^}JVo) = ^T^ / - / exp{-/3[F// 



i{Xi}N) 



Vrnf{{Xt}N;{yz}No)- / ^J'ix) p{x)]} , 



(21) 

where V//({a;i}jv) and Kn/dxjjv; {yjwo) are the to- 
tal fluid-fluid interaction potential and the matrix-fluid 
interaction potential, respectively, and p{x) is the mi- 
croscopic density of fluid particles at point x. The 
Helmholtz free energy of the fluid (or effective action in 
fleld-theoretical language) is as usual the Legendre trans- 
form of the grand potential flf{p,). 

To perform the average over the matrix configurations 
of the logarithm in Eq. d^U]) . one can use the replica 
trick^ which simply states that 



Bnf(n) = lim - 



exp[nlogS/(^;{j/,}7Vo] + l) 

li-^(l-exp[-/3fi7^(^;7Vo)]) (22) 

(3nyPip;Noy 



lim 

n-s-O 



where 



equations"—: 



/3fi7^(M;iVo) = log 



E 



1 



^'"o Ni,...,Nr 



Ni\...NJ. 



/ - / ^o({y»}A'o)exp{ -/3^ [V7/({a;„}A'„ 

+ Vmfi{Xaz}N^]{yx}No) - / tJ'{x)Pa{x)]} , 



with Pa{x) the microscopic fluid density of rephca a. 
rff^{p;NQ) corresponds to an equilibrium thermody- 
namic potential for a mixture of n fluid replicas at the 
chemical potential fi and one matrix species at fixed num- 
ber of particle iVo (and it is proportional to the number 
n of replicas) . Thanks to the replica trick, the quenched- 
annealed binary mixture has been replaced by a fully 
annealed mixture of n -I- 1 species. The fluid particles do 
not interact directly when they belong to distinct repli- 
cas but they all interact with the matrix particles. In 
the following, the (now annealed) matrix is character- 
ized by a subscript and the fluid replicas take the in- 
dices l,2,...,n. We denote by Greek letters, a,/3, ..., the 
n + 1 species {e.g., a = 0, 1, ..., n) and by Roman letters, 
a, b, ..., the n fluid replicas {e.g., a ~ 1, ..., n). 

One can then apply all the tools of liquid-state the- 
ory to this (n -|- l)-component equilibrium mixture. In 
particular, one can introduce direct and total correla- 
tion functions (or in field-theoretical language, proper 
or IPI vertices and Green's functions) and relate them 
via OZ equations that stem from the Legendre trans- 
form between grand potential (generating functional of 
the Green's functions) and Helmholtz free energy (effec- 
tive action, which is the generating functional of the IPI 
vertices). At the pair level and using the field-theoretical 
notations (see section [ll]), one has in Fourier space and 
in matrix form 



W^^\q)=T(^^\q)-^ 



or equivalently in liquid-state notations. 



p + ph{q)p= [p ^ -c{q)] 



(24) 



(25) 



where the matrices are (n + 1) x (n + 1) and p is diagonal 
with Poo = Po and pn = ... = /?„„ = p. 

To recover the properties of the original system, an an- 
alytic continuation to arbitrary values of n together with 
the limit n — > have to be taken. Assuming that the per- 
mutational symmetry among the n replicas is not broken 
(contrary to what occurs for instance for the Sherrington- 
Kirkpatrick spin-glass model^), one can easily diagonal- 
ize the (n + 1) x (n + 1) matrices and one ends up in the 
limit n — >■ with the following "replica-symmetric OZ 



hmf{q) = Smrn{q)- — T^ 

1 - pCcon{q) 



1 + phcon{q) = 

hdis{q) 



1 - PCcon{q) ' 
Cdis{q) + PoSm7n{q)Cmf{q)'^ 
[1 - pCcon{qW 



(26) 



In the above equation, Smm{q) is the structure factor 
of the matrix phase of the porous material, /icon and 
hdis are the so-called "connected" and "disconnected" 
(or "blocking"™) fluid-fluid pair correlation functions. In 
physical terms, the latter are defined from the 1-body 
and 2-body densities Pf{x) and pff{x,y) of the fluid in 
a given configuration of the solid matrix as 



P hcon{\xi -X2\) = Pff {xi,X2) - Pf{xi)pf{x2) 
P^hdis{\xi -~X2\)^ Pf{xi)pf{x2) ~ P^, 



(27) 



and their sum is equal to the total fiuid-fiuid pair correla- 
tion function hff that is directly related to the fiuid-fiuid 
structure factor through Sff{q) = 1 + phff{q); finally, 
p — 'pj is the mean fiuid density. 

In terms of the replicated mixture, one has 



hcon{q) = lim haa{q) " hdis{q) 

n— >0 

hdis{q) = lim hab{q), 

rn-O 



(28) 



where a ^ b in the last equation. Similarly, the 
static structure factor of the matrix can be expressed 
as Smm{q) = 1 + po lim„^o /ioo('z) and the matrix- 
fiuid correlation function as hmf{q) = lim„_>o /^Oa (?) — 
\iuin^ohao{q)- 

The connected fiuid-fiuid pair correlation functions and 
the average Helmholtz free-energy density for the fiuid in 
the porous material Af{p;pQ) are related by a compress- 
ibility sum rule^. 



Cco7i{q = 0:p) = 

P 



1 



p + p^hcon{q = 0;p) 
^ d'l3Af{p;Po) ^''^ 

A similar relation holds for the matrix-fiuid direct corre- 
lation function: 



'Cmf{q = 0;p) = 



d^PAf{p;po) 
dpdpo 



(30) 



However, neither the total nor the direct disconnected 
fiuid-fiuid correlation functions are related to derivatives 
of the average Helmholtz free-energy density for the fiuid 

Af{p). 

Over the last two decades, there has been a great 
deal of activity to solve the replica-symmetric OZ equa- 
tions with various standard closures'*^ (Percus-Yevick, 



HNC, MSA, ORPA, EXP and OCT, etc..) and vari- 
ous choices of microstructure for the sohd matrix [which 
enters the equations only through the structure fac- 
tor Smmiq)]- Fluid phase diagrams and pair correla- 
tion functions have been computed and compare well, 
at least semi-quantitatively, with the existing simulation 
and experimental data^s^"— . However, the standard 
approximations break down near the critical points, when 
present, because of their inability to describe correctly 
the long-range fluctuations. A promising route is then to 
apply the HRT to the quenched-annealed mixtures. 



IV. A NAIVE HRT APPROACH 



leads to 



A. Strategy 



dtrrHp.}] = 



¥' 



tri log 



-,rep(2) 



[{Pc}] 



\{q,-q), 



(33) 



where tr is a trace over the n -I- 1 components. The 
whole HRT hierarchy for IPI vertices (direct correlation 
functions) evaluated for uniform density fields is again 
simply obtained from the above functional equation by 
repeated differentiations. 

To put this formal scheme at work, we have to (i) con- 
sider an analytic continuation in the number n of replicas 
and take the limit n — > 0, (ii) devise an approximation 
to truncate the HRT hierarchy, and (iii) choose a specific 
cutoff function Rk{q)- 



The strategy that appears natural to describe the crit- 
ical behavior of fluids in a disordered porous medium 
consists in generalizing the HRT approach summarized 
in section |ll] to the replicated mixture with n + 1 com- 
ponents introduced in section IHII We choose for ref- 
erence system the mixture defined with all interaction 
pair potentials (matrix- fluid and fluid-fluid), except the 
attractive interaction Wff between fluid particles in the 
same replica (recall that there are no direct interactions 
between fluid particles belonging to different replicas). 
Generalizing Eq. ^, we therefore introduce a modified 
interaction matrix Wk-ap in which the only nonzero el- 
ements are the diagonal terms Wk-aa with a = l,...,n, 
which are defined such that 



Wk;aa{q) = Rk{q)wff{q), 



(31) 



with Rk{q) having the same properties as for a one- 
component fluid (see section [H]) . As before, we also in- 
troduce (j)k{q) = -(iRk{q)wff{q). 

With the above choice of reference system and IR regu- 
lator, the extension of the exact functional HRT equation 
to the equilibrium {n + l)-component mixture reads 



dtTriip.}] 

2 / dtMq) 



E 

a=l 



[Tr''\{p.}]+4>ul 



(9. -q) 

(32) 

where we have momentarily introduced inhomogeneous 
density fields Pa{x) that are different for the different 



replicas; ^ ^ 

k-afi 
(^0/3(1 - 5ao) 



rep(2) 



and 



with rz^':>{x,y) 



1 are (n + 1) x (n + 1) matrices 
S'^Tl''^ /Spa{x)6pi3{y) and lap = 
As stressed already several times. 



the 



effective average action F^ '^ is, up to a factor /3, the 
Helmholtz free-energy functional and the 2-point IPI ver- 
tices F^^^l are related to the direct pair correlation func- 
tions [see Eqs. (|^ and d^ ]. Another more compact way 
to write the HRT equation is by introducing an opera- 
tor dt that acts only on the fc-dependence of 4>k {q) ■ This 



B. Replica-symmetric HRT 

Concerning the first point mentioned above [point (i)], 
we make the assumption, which is common in the treat- 
ment of quenched-annealed mixtures (see section IHip . 
that the permutational symmetry between replicas is not 
broken when n — > (provided all the densities pa are 
taken as equal, which they should be when describing 
the physical system). This implies in particular that any 
{n + 1) X {n + 1) matrix M appearing in the HRT hi- 
erarchy has the same structure, with Afoa = -^^01 f^nd 
Mab = (Mil - Mi2)Sab + M12, for a,b = l,...,n. Such 
a matrix can be easily diagonalized, which henceforth 
allows for straightforward algebraic manipulations. For 
instance, one simply obtains 



log(A/n - M12) 



M12 



Moo 



M- 



11 



M- 



12 



<rlogM:=log(A/oo)H 

+ 0{n^). 

(34) 
Similar but more complicated properties hold for the 

higher-order tensors, e.g., Tl':^f^{x,y, z). 

The first equation of the HRT hierarchy in the n — s> 
limit then provides a flow equation for the (modifled) 
Helmholtz free-energy density of the fluid at scale k, 
-4/,fe(p;po) = \imn~^od[Tl''^{p;po)/V]/dn, which reads 



dtAf^kip; Po) =7: / dt(j)k{q) 



1 



^con.k yq ) 
Cdis,k{q) + PoSmm{q)Cmf,k{qY 



[P-' 



k{qW 



(35) 
where the direct correlation functions at the scale k de- 
pend on p and po a-nd are related to the IPI vertices 
according to 

p-' - c^ouAi) = ri:ff \g) - Tlj!^\q) + Mq), 

Q..,.(g) = -F;jf(g), (36) 

Cmf.kil) = -rfc;of' (9) 



in the limit n —^ 0. The first term of the right-hand side 
of Eq. dSS]) is similar to that in the HRT for the bulk fluid 
[see Eq. ^^] and the second term is really the signature 
of the presence of the disordered porous matrix: it van- 
ishes identically when po = while, then, Ccon.k = Ck- 
Note that the total pair correlation functions are related 
to the above direct correlation functions through an ob- 
vious generalization at scale k of the replica-symmetric 
OZ equations, Eq. ((26)) . 

The other equations of the hierarchy can be written 
in the same way as Eq. psp . The derivation is quite 
straightforward, but, because of the rapidly increasing 
number of irreducible components of the many-body cor- 
relation functions that one must consider when deal- 
ing with higher orders, the explicit expressions are too 
lengthy to be reproduced here, even in a diagrammatic 
representation. (An illustration will be given below for 
the flow equation of the connected pair correlation func- 
tion in the asymptotic regime near the critical point.) 
Note that there are no flow equations for the pure ma- 
trix functions in the limit n ^> 0, which is due to the fact 
that the matrix is rigid and independent of the fluid. 

Exact compressibility-like sum rules are also obtained, 
that relate the connected component of the fluid-fluid di- 
rect correlation function and the matrix-fluid direct cor- 
relation function at zero wavevector to the second deriva- 
tives of ^/,fe: 



1 

P 



,k{q = 0) - Mq = 0) 



d^Af,k{p;po) 
5p2 



-Cmf,k{q = 0) 



d^Af^k{p\Po) 
dpdpo 



(37) 



(38) 



Generalizations of these identities can be derived for the 
higher orders. 

C. Asymptotic behavior 

When approaching the gas-liquid critical point (in the 
cases of course where the effect exerted by the porous 
material is not so strong as to destroy the gas-liquid 
transition of the adsorbed fluid), the fluid compressibil- 
ity diverges. Accordingly, the total connected fluid-fluid 
correlation function at the critical point behaves as 



'^con,crit\Q ) ^^ Q 



(39) 



at small wavevector q. The asymptotic regime of the 
HRT equations describes the long-distance physics near 
the critical point, which corresponds to small IR cutoff 
k and small wavevectors. Keeping in mind the definition 
of the total and direct correlation functions, the above 
critical behavior translates into the following small k be- 
havior of the connected pair correlation functions: 



1 



- Ccon,k{q) - ^fc(<?) — fc^ '' Mcon,fe(g^), (40) 



where q^ = q^ /k"^ and the scaling function Ucon,k is ex- 
pected to converge to a finite limit Ucon,* when A; — > 0. 
To be consistent with this behavior, we choose the IR 
regulator such that 4>k{q) — k'^^^; more specifically, we 
define 



q 



Rk{q) = Zkk'ri^), 



(41) 



where the function r(y) is equal to 1 when y = and 
decreases fast enough to zero when y — ;■ cx) [to cite 
a few examples used in the nonperturbative-RG liter- 
ature: (1 - y)H{l - 2/)49, y/(exp(2/) - 1)29,30^ ^^ ^^^ 
sharp cutoff H{1 — y), where as before H is the Heav- 
iside step function]. Z^ ~ k~^ is introduced to account 
for the so-called field renormalization (which introduces 
the "anomalous dimension" 77) and is defined through a 
prescription on the direct pair correlation function (see 
below). In the asymptotic regime, for wavevectors such 
that q/k — q = 0(1), the cutoff interaction 4>k{q) can 
then be replaced by 



M<l)--l3wffiq^0)Zkk\{f), 



(42) 



where we recall that Wff(0) is negative. One then finds, 
after dropping some irrelevant constants that can al- 
ways be incorporated in a redefinition of the functions 



and variables, that p 
-2 



- Ccon.kil) ~ fc ''[Ucon,fe(g ) + 

r{q'^)], which gives back Eq. ([55)1 (provided again that 
Ihnk-^o Ucon.kiq'^) is finite and different from zero). 

On the other hand, at the critical point, the matrix- 
fluid direct correlation function at zero wavevector satis- 
fles 



Cmf,crit{q = 0) 



d'Af 



dpdpo 



dpc 



dpo 



(43) 



which is finite. Therefore, in the asymptotic regime, 
PoSmm{q)cmf,k{q)^ is finite. The fiuid-fluid discon- 
nected total pair correlation function hdis.k{q), which is 
given in Eq. (|26p . has then two possible types of be- 
havior. If an accidental cancellation occurs such that 
Cdis.k{q) + PoSmmiq)cmf,k{q)'^ -^ 0, then the two terms 
of the right-hand side of Eq. ([55)1 could possibly scale 
in the same way if Cdis.k{q) + PoSm7n{q)cmf,k{qY goes 
to zero as fc^-i)^ This is actually the case when a 
special symmetry is present that leads to a cancella- 
tion of the random chemical potentials exerted by the 
solid matrix on the fluid. For a lattice-gas version of 
the quenched-annealed mixture, we showed that when 
the ratio between the attractive matrix-fluid interac- 
tion and the fluid-fluid one is exactly equal to 2, the 
system is isomorphic to the so-called site-diluted Ising 
model^°. In such a case, the symmetries of the Hamilto- 
nian impose that, for the fully interacting system (i.e., 
k = Q) along the critical isochore, Cdis,k=o{q) = whereas 
Cmf,k=o{q) = (1/2) [-(0"^ -I- Ccon,k=o{q)] gocs to zero as 
one approaches the critical point^. As a consequence, 

Cdis.kil) + PoSmrniq)Cmf,k{l)^ -> aS l-pCcon,kil) ^heu 

A: ^' in the critical region. The site-diluted Ising model 



has an upper critical dimension of 4 and quenched dis- 
order plays a much weaker role there than in the RFIM. 
It should thus be handled in a specific manner (in par- 
ticular, if one wishes to preserve the symmetry property 
along the flow, one needs to make the IR regulator also 
act on the matrix-fluid attractive interaction) and will 
not be further considered in this paper. 

In the generic case, one expects Cdis,k(<l) + 
PoSmm{q)cmf ,k{<lY ¥" 0- The second term of the right- 
hand side of Eq. (|35|) then asymptotically dominates the 
flrst term at small fc. As a result, the latter is irrelevant 
and can be dropped when fc — ^ 0. The quenched disor- 
der then controls the critical behavior. Technically, this 
corresponds to a "zero-temperature fixed point" which is 
characterized by an additional exponent 9 and a violation 
of the so-called hyperscaling relation. 



2-a = v{d~e), 



(44) 



where a and v are the specific-heat and correlation-length 
critical exponents. Such a zero-temperature fixed point is 
known to control the critical behavior of the RFIM^"— . 
One then anticipates that the disconnected component 
asymptotically behaves as 



Cdis,k{q) 



-(277-77) 



Udis,k{q^ 



(45) 



with the exponent i) related to the exponent 9 through 
f] = 2 -\- T] — 9 and satisfying ij < fj < 2r]. 

In the specific case where the upper bound is satu- 
rated, f] = 2?7, the disconnected direct correlation func- 
tion goes to a constant, just as the matrix-fluid contribu- 
tion and the two terms must be kept in the expression of 
hdis,k{q)- On the other hand, if ^ < 277, the matrix-fluid 
contribution is subdominant when compared to the dis- 
connected fluid-fluid one in the asymptotic regime and 
can be dropped. This is also true at the higher orders of 
the HRT hierarchy. 

Note that the emergence of a second "anomalous di- 
mension" 77 leads to the following behavior of the discon- 
nected pair correlations at criticality: 



hdis,crit\<lj ^ Q 



-{i-n) 



(46) 



which implies that the static structure factor of the fluid 
Sffiq) = ^ + phcon{q) + phdzsiq) also diverges as g"*"'"''), 
i.e. much more strongly than in the bulk fluid. 

To express the HRT hierarchy in the asymptotic 
regime, it is convenient to introduce "dimensionless" 
quantities. Indeed, it is only by eliminating the scal- 
ing dimensions that the scale-free behavior associated 
with the existence of a fixed point can be found in the 
RG flow. Near a zero-temperature flxed point, the con- 
ventional scaling dimensions are modificd^^ due to the 
presence of the new exponent 9 (or equivalently fj), and 
one introduces ip — fc~'''~^"'"''"^(p — Pcrit), CLki'p) = 
/j-(rf-«)/2(_4^^^(p) _ Af.k-crit) (where the subscript crit 
indicates that a quantity is evaluated at the critical 
point), etc. All constants, such as —(3wff{q = 0) in 



Eq. p2|) , can be included in a redefinition of the field as 
well as of the various functions and vertices; they will no 
longer appear explicitly. 

The asymptotic form of the replica-symmetric HRT hi- 
erarchy is derived with the help of above considerations 
and by dropping all subdominant terms. The first equa- 
tion then reads 

dtak{(p) + (d- 2 - T] + fi)akiip) - -{d - 4: + f])ipd^ak{Lp) 



^J[{2-r^)rif)-2cfr\f)]- 



udisAq'^^v) 



[Wcon,fe(9^;<^) +r(^2)]2 

(47) 

where r'{y) is the derivative of r(y). The above equa- 
tion is valid if ^ < 2?]. In the case where fj = 2r], one 
must also take into account the matrix-fluid contribution, 
PO'S'm?n(0)c^i Tj, that is no longer irrelevant compared to 
Cdis,k- An illustration of the higher-order equations of 
the HRT hierarchy is given in appendix \K\ 



D. Replica-symmetric HRT for the RFIM 

The random-field Ising model (RFIM) on a d- 
dimensional lattice is described by the following 
HamiltonianSii^: 



H — —J 2_^ SiSj — 2_^ hiSi, 



<ij> 



(48) 



where Si = ±1 are Ising variables, < ij > indicate a sum 
over distinct pairs of nearest neighbors, and the random 
field hi is delta-correlated in space and is drawn from a 
symmetric probability distribution with zero mean and 
variance As, i.e., hi — 0, hihj — SijAs. Choosing 
for simplicity a Gaussian distribution and applying the 
replica trick, one finds that the properties of the model 
can be obtained from the "replicated" Hamiltonian 



jjr 



Af 



^-JJ2Y1 '^"•5'ja - 1^ Yl XI '^^°'^*'" (49) 



a=l <ij> 



a, 6—1 i 



with, as before, the n ^- limit taken at the end of 
the calculations^^. Note that the replicated model has a 
global Z2 symmetry, the Hamiltonian being invariant by 
the simultaneous change Sia -^ —Sia in all replicas. This 
symmetry is absent in the quenched-annealed mixture at 
the microscopic level. 

One can repeat the steps of the replica-symmetric HRT 
approach with a change of notations from fiuid to mag- 
netic language: the density is replaced by the magne- 
tization m=< Sia >, the Helmholtz free-energy density 
by the Gibbs free-energy per spin Akirn), etc. The 
reference system is that with only on-site interactions, 
-{A.B/2T)J2'2,b=iJ2iSiaSib, and the Ising constraint 
^ia — 1 (which is equivalent in a lattice-gas language 
to the "core" or single-site occupancy condition). The 
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attractive potential is simply —JXij where A^ is the ad- 
jacency matrix of the lattice so that its Fourier trans- 
form X{q) — J2u=i cos(g,/) for a d-dimensional cubic 
lattice. An IR regulator is introduced and one defines 

The HRT hierarchy is then formally the same as that 
derived above for the quenched-annealed mixture, except 
that all the direct correlation functions involving the ma- 
trix disappear from the flow equations. For instance, the 
first equation of the HRT hierarchy is simply 



dtAk{m) 



1 



2 / dt4>k{<i) 



,k{q;m) 



Cdis.ki'l','^) 

Ccon,k{q]rnY 



(50) 



where the direct correlation functions are conventionally 
defined in a slightly different way as in liquid-state sta- 

(2) 

tistical mechanics with Ccon,k{<l', m) = F^ {q; m) + (j)k{q) 
(on the other hand, Cdis,k{Q] m) has the same definition). 
With the correspondence Ccon,k{q; m) -^ 1/ p-Ccon,k{.q] p) 
[and Cdis,k{q\m) -^ Cdis,k{q] p) + pSmm[q)cmf,k{q\ p)\ 
Eq. ([50| coincides with Eq. (|35|) . This carries over to the 
higher orders. 

The critical behavior of the RFIM is known to be 
controlled by a zero-temperature fixed point^^— and 
the same scaling as in section IIVCI is used to put the 
flow equations in a dimensionless asymptotic form (of 
course, there is no need here to consider the matrix- 
fluid functions) . These asymptotic flow equations for the 
RFIM are then identical to those derived above for the 
quenched-annealed mixture. (They are also identical to 
the flow equations obtained in the field-theoretical for- 
mulation of the RFIM.) 

To make fully rigorous the proof that the critical be- 
havior of the quenched-annealed system and that of the 
RFIM are identical, one must ensure that the initial con- 
ditions at the microscopic scale of the two systems do 
lead to the same fixed point. As previously noted, the 
symmetry of the initial condition of the RG fiow is dif- 
ferent in the two cases since the fiuid does not have the 
(statistical) Z2 global inversion symmetry of the Ising 
model. However, provided that the asymptotic HRT 
equations (which are formally the same for the RFIM and 
the quenched-annealed mixture) admit a unique fixed 
point with the appropriate properties to describe a crit- 
ical point, the fiow equations should lose the memory of 
the initial conditions and the critical behavior of the two 
systems be the same. (The same is true for the critical 
behavior of bulk fluids which is in the same universality 
class as that of the pure Ising model^ despite the absence 
of Z2 symmetry.) The fact that the (critical) fixed point 
is unique will be illustrated with specific approximations 
below. It is also worth pointing out that the above rea- 
soning, as previously discussed, assumes that no acci- 
dental cancellation takes place in the quenched-annealed 
case between disconnected and matrix-fiuid components, 
which should also be checked in actual computations. 



E. Ornstein-Zernike approximation 

Closures of the HRT hierarchy for the quenched- 
annealed mixture can be devised by extending the ap- 
proximation schemes developed for the HRT of bulk flu- 
ids (see section|TT| to the replicated equilibrium mixture. 
The simplest closure is an RPA-like approximation (see 
section |ll| that enforces the compressibility sum rules, 
Eqs. (PI and (PI): 



— Cco„,fc(g; p) - (pk{q) = V^ 

p op"^ 



Cmf,k{q;p) = - 



d^Af,kip;po) 
dpdpo 

+ [c,nf,Riq;p) 



•»/,«(? = 0;p)]' 



Cdis,k{q; p) = Cdis,R{q = 0; p). 



(51) 



(52) 



(53) 



where as before the subscript R indicates a quantity cal- 
culated in the reference system in which the fluid atoms 
interact with the matrix with the full pair potential and 
interact among themselves through the steep repulsive 
interaction component. [Eq. ([51]) for the connected func- 
tion could also be replaced by an equation similar to 
Eq. ([HI).] 

As for the bulk fiuid, the above closure could be im- 
proved to better account for the physics at short dis- 
tances. For instance, if the steep repulsive pair interac- 
tions between fluid atoms and between fluid and solid 
atoms are modeled as hard-core potentials, one can use 
an ORPA-like closure by adding contributions to the di- 
rect correlation functions for distances less than the rele- 
vant core diameter: only two core conditions however ex- 



ist, /i//,fe 



''dis.k 



-1 inside the fluid-fluid core 



range and ft//.fc = — 1 inside the fluid-matrix core range. 
Note also that, due to the absence of compressibility sum 
rule relating the disconnected correlation function to the 
fluid Helmholtz free-energy density, there is no obvious 
way to improve the closure for this correlation function 
beyond Eq. ([53]) . As a result the disconnected direct cor- 
relation function is not renormalized in any of the above 
approximations . 

Irrespective of the degree of refinement, most common 
approximations of the direct correlation functions are of 
OZ type, which means that their small-wavevector be- 
havior is regular, i.e., 



- - ccon,k{q\ p) - (t>k{q) = — -^-2 ^ 0{q ) 

p Op'^ 



Cnif,kiq]p) 



d^Af,k{p;po) 



O(g'), 



(54) 



dpdpa 
Cdis.kiq; p) = Cdis,R{q = 0;p) + 0{q^). 



From the regular g-dependence of these direct correla- 
tion functions, one immediately derives that the critical 
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exponents 77 and rj satisfy 

r] = fj = 0. (55) 

Generically with such an OZ approximation, 

p^^ - Ccon,k{q; p) -4>k{q) -d'^Af^k{p;po)/dp'^ - buq^ 

(56) 
and 



Cdis,k{q) + PQSmm{q)Cmf ,k{qY — ^k 



(57) 



in the small-fc, smah-g limit, where 6^ > and K^ > 
are regular functions of fc, p, and T (and of the disorder 
characteristics such as po) with finite limits, 6 > and 
5 > when fc — >■ 0. (Note that as ry = 77 = 0, the matrix- 
fluid contribution is not subdominant with respect to the 
disconnected one.) 

The asymptotic HRT equation, Eq. (|T7l). can be reex- 
pressed as 



dtak{ip) + (rf- 



2)afc((p)--(rf-4)(^4((p) 



2[r{f)-q^r'{f)] 



(58) 



where a prime always denotes a derivative with respect 
to the argument and the constants b and K have been 
incorporated into a trivial redefinition of a^ and ip. It is 
easy to show that the very same equation as Eq. ([55)) is 
obtained under the same assumptions for the RFIM. In 
the context of field theory, this kind of approximation to 
the RG fiow is known as "local potential approximation" 
(LPA)^i^. 

For actual computations, we have considered the sharp 
IR cutoff (see section |lT| , so that the asymptotic HRT 
equation becomes 



dtak{^) = - (rf - 2)ak{ip) + i^{d - 4)^0'^ (93) 



1 + <(¥>) l + <(^ = 0) 



(59) 



The fixed point that describes the gas-liquid critical 
behavior of the fluid is determined by looking for the 
scale- free (stationary) solution of the above equation, i.e. 
dtak\* = 0. Taking for convenience one derivative with 
respect to the field, we then obtain 



d<(^)-(d-4)^<'(^) 



<'{V) 



[1 + <M]' 



0. 



(60) 



One has to look for a solution that behaves at large pos- 
itive and negative field (</? — > ibcx)) as a'^((p) '^ |(^|'^/('^~'*) 
and that is defined and continuous over the whole real 
axis. This condition ensures that the second-order dif- 
ferential equation has only a small number of acceptable 
solutions^. 

When the spatial dimension d > 6, it is easily shown 
that the only fixed-point solution is a^(ip) — 0, which 



corresponds to the Gaussian fixed point. The exponents 
are then the classical ones obtained from mean-field the- 
ory; d = 6 is therefore the upper critical dimension of the 
system at criticality, which represents a shift by 2 com- 
pared to bulk fluids and results from the strong effect 
of the quenched disorder (random chemical potential or 
random field). Below 6, a new nontrivial solution with 
properties associated to a critical point is found in addi- 
tion to the Gaussian fixed point, and the latter becomes 
unstable. 

The fact that the anomalous dimensions rj and 77 
are zero in the present approximation prevents us from 
studying dimensions below 4: as can be seen from 
Eq. ([^5)1 . which is equivalent in real space to hdis.critir) ^ 
j,-(d-4+fj)^ the critical decay of the pair correlations is 
not ensured in this case and, in fact, no critical point 
is expected. We have therefore calculated the solutions 
of Eq. ([60]) in d = 5 by using a shooting method as in 
ref. [57] . We have found a unique nontrivial solution 
a'^{(p) satisfying the conditions described above, and it 
is an odd function of the field. The corresponding a*((/3) 
is then an even function. As for the bulk fluid whose 
critical behavior is in the same universality class as the 
Ising model, one recovers that the Z2 inversion symmetry 
(which for liquids amounts to a symmetry with respect to 
the critical density) is obeyed by the critical theory, as in 
the RFIM, despite the fact that the initial Hamiltonian 
of the quenched-annealed model is not symmetric. 

The stability of the fixed point and the critical expo- 
nents (if the fixed point indeed properly describes the 
critical behavior) are investigated by linearizing the fiow 
equation around the fixed-point solution and studying 
the spectrum of eigenvalues. Because of the asymptotic 
Z2 symmetry of the fixed point, the eigenvalues can be 
sorted according to the odd or even behavior under field 
(fi reversal of the associated eigenfunctions. The largest 
eigenvalue Ae corresponding to an even cigcnfunction de- 
scribes the relevant direction associated with the temper- 
ature or the disorder strength whereas the largest eigen- 
value Ao associated with an odd eigenfunction describes 
the chemical-potential (or magnetic-field for the RFIM) 
dependence along the critical isotherm. All other eigen- 
values should be negative for a proper critical fixed point. 
Writing ak{f) — a*((/9) -|- k'~^ff^[ip) and then linearizing 
around the fixed point lead to 



(-A + d-2)/A(^)-^¥'/;M 
/a(0) 



m^) 



[1 + <(^)]2 



[l + <(0)] = 



(61) 

The "odd" eigenvalue Aq is obtained analytically by 
noticing that the above equation has an odd solution 
/a((/9) = a'^[(p) with Ao — (d ~ 4)/2 [compare with 
Eq. ([60t ]. The "even" eigenvalue Ag and its associ- 
ated eigenfunction on the other hand must be numer- 
ically determined. (One looks for an even solution of 
the second-order differential equation that is, as before. 



12 



defined and continuous everywhere but behaves now as 

^|2(d-2+A)/(d-4) ^j^gj^ ^ _^ QQ ■) One then obtains 

from these eigenvalues the following critical exponents 
in d = 5: 



v= — ^ 0.6496, 

A, 



7 
5 - 



2v ^ 1.299, 

2 _ d 
An^ d- 



(62) 



1 



where 7 is the exponent describing the divergence of 
the connected susceptibility {d^Af/dp^)~^, whereas 5 
characterizes the relation between density and chemi- 
cal potential along the critical isotherm. (Recall that 
77 = fy = so that 9 = 2; one can also define an exponent 
for the divergence of the "disconnected susceptibility", 
7 = (4 — ^)z/ ~ 2.598.) In addition, we have computed 
the second largest eigenvalue corresponding to an even 
eigenfunction. It is negative, as required, and its abso- 
lute value provides the critical exponent uj controlling the 
corrections to scaling. The result is 



0.6557. 



(63) 



It is puzzling to find that the exponents predicted by 
the above OZ approximation for the critical behavior of a 
fluid in a disordered porous matrix (and equivalently the 
LPA for the RFIM) do not obey the so-called "dimen- 
sional reduction"^. According to the latter property, the 
critical behavior of the RFIM in dimension d is the same 
as the critical behavior of the pure Ising model without 
quenched disorder in dimension d — 2. The critical expo- 
nents for the system in the absence of quenched disorder 
(bulk fiuid or pure Ising model) in d = 3 are equal to 



0.689,7-1.378,(5 = 5, 



0.581 



(64) 



in the OZ/LPA approximation with a sharp IR cutofF^'^. 
Except for S, they are thus different from the exponents 
obtained above for d = 5 (even when considering error 
bars). 

The dimensional-reduction property of the RFIM is 
predicted by perturbation theory at all orders^"—, but 
has been rigorously proven wrong in low enough di- 
mensions {e.g., in d = 3) ^^'^'^ . However, from pre- 
vious studies on the RFIM, one does not expect that 
the breakdown of dimensional reduction could be cap- 
tured by the simple replica-symmetric formulation used 
above^^i^. One indeed finds that dimensional reduction 
is verified at first order in the expansion in e = 6 — d 
around the upper critical dimension in the OZ/LPA ap- 
proximation: in this case, one easily obtain from Eq. (|60|) 
that a^,{(p) oc e{(fi'^/2 — ip'^/12); linearizing around this 
fixed point leads to 



1 
"=2 



12'^ 



1 



h^ 



(65) 



which exactly corresponds to the result at first order 
in the expansion in e = 4 — d around the upper crit- 
ical dimension of the pure Ising model or the bulk 



fiuid^. At second order in the e-expansion however, 
the dimensional-reduction property is lost in the present 
approximation^ while it is still verified in the exact per- 
turbative expansion. Clearly the breakdown of dimen- 
sional reduction in the present approach is a consequence 
of the approximations made, not a genuine prediction. 



V. FLAWS OF THE NAIVE HRT AND 
LESSONS FROM THE RFIM 

A. Flaws of the naive HRT approach 

The naive HRT approach presented above has three 
major flaws: 

(1) The fact that dimensional reduction is not obeyed, 
despite the simple nature of the approximation that 
cannot describe the complex physics associated with 
dimensional-reduction breakdown (metastable states for 
instance^), is likely to be an artifact of the chosen ap- 
proximation. In the case of the RFIM, dimensional re- 
duction has been shown to follow from an underlying 
"supersymmetry" of the field theory describing the long 
wavelength behavior at criticality^^. As we have obtained 
indications that the critical behavior of a fluid in a dis- 
ordered porous material and that of the RFIM are the 
same, we can anticipate that the same underlying super- 
symmetric property could in principle be asymptotically 
found in the quenched- annealed mixture. Clearly, one 
should then distinguish between a spontaneous break- 
ing of the supersymmetry, which is the physical mech- 
anism for the (necessary) breakdown of dimensional re- 
duction in low dimensions (at least in d = 2, 3), and an 
explicit breaking of supersymmetry due to an inappro- 
priate choice of approximation and cutoff functions. The 
naive approach most likely breaks the supersymmetry ex- 
plicitly. However there is no way to detect, and further to 
cure, this problem within the formalism so far developed. 

(2) The OZ/LPA approximation does not allow for 
fleld renormalization, which implies that rj = fj = 0. 
As the spatial dependence of the 2-point Green's func- 
tion/total pair correlation function at the critical point 
are predicted to go as r^'''~*+^) for the disconnected 
component, no criticality is found for dimensions d < 4 
and one cannot directly study the physically relevant 
three-dimensional situation. It is known that going be- 
yond the OZ/LPA approximation with a hard-cutoff for- 
mulation is difficult due to the proliferation of singulari- 
ties associated with derivatives of the cutoff function. 

(3) The replica-symmetric formalism, besides being 
wrong when a hypothetical spontaneous breaking of the 
permutational symmetry between replicas takes place, 
does not allow one to easily close the equation at the 
level of the disconnected 2-point vertices. As we have 
seen, there is no compressibility sum rule relating the 
latter to the average Helmholtz free-energy functional of 
the fluid. As a result, we had to assume that the discon- 
nected fluid-fluid direct pair correlation function is not 
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renormalized and essentially remains equal to its values 
in the reference system. (On the other hand, the matrix- 
fluid direct pair correlation function can be renormal- 
ized but the effect of this renormalization is benign as 
d'^Af /dpdpQ stays finite at criticality.) 

Going beyond the naive HRT thus requires to find ways 
for addressing these shortcomings. 



tween IPI vertices (direct correlation functions) that are 
known as Ward-Takahashi identities'^^. For instance su- 
persymmetry in the RFIM implies that asymptotically, 
i.e. near the critical point and at small momenta (long 
wavelengths)^", 



.(11) 

fe2 



{q^]m,m) = l^dq2T)^{{q^;m) 



(66) 



B. A detour via the RFIM and superfield theory 

It actually took a long detour via the (super)field 
theory of the RFIM and its nonperturbative RG 
description22r^ to make progress. To keep a long story 
short we only briefly summarize the lessons that can be 
drawn from the latter study. 

First, the replica-symmetric formalism should be re- 
placed by a framework in which the replica symme- 
try is explicitly broken by letting the sources (in field- 
theoretical language) or chemical potentials (in fluid lan- 
guage) be different for each replica. This a priori unphys- 
ical situation (the fluid in the disordered porous medium 
or the Ising model in a random field are subject to a 
unique chemical potential or external magnetic field) is 
actually a means to properly generate the "cumulants of 
the renormalized disorder" . Indeed, in the presence of 
quenched disorder (and before any use of a replica trick), 
the generating functionals and thermodynamic potentials 
of the fluid (or Ising model) are random quantities. As 
such, they should be described by a functional probability 
distribution or by an infinite set of cumulants, the lat- 
ter being more convenient as the cumulants are averaged 
quantities for which statistical invariance under trans- 
lations and rotations is recovered. It is easily realized 
that the only way to describe the cumulants of order 2 or 
more with their full functional dependence is by consider- 
ing copies or replicas of the initial system with the same 
disorder but different density fields, hence different chem- 
ical potentials (or sources) . This point is central as it has 
been found in the case of the RFIM that the influence of 
the rare collective events stemming from the presence of 
quenched disorder (abrupt changes of the ground-state 
configuration under a variation of the external source 
known as "avalanches" or "shocks" and near-degenerate 
metastable states known as "droplets") precisely results 
in a singular behavior of the second and higher-order cu- 
mulants in their functional dependence^^Ti^i. This will be 
made more explicit below. 

Secondly, the underlying supersymmetry does not 
show up easily in the conventional field formalism, even 
when one considers an explicit breaking of the replica 
symmetry. It is made explicit when one upgrades the 
description to a superfield approach that arises for for- 
mal manipulations of the original problem''^. There is 
no point to go here through the derivation. It suffices 
to say that the supersymmetry can be viewed as an 
asymptotic property of the theory and that its conse- 
quences in the present formalism are exact relations be- 



with the obvious notation: 



.(11) 



{2T:Y5^'^\qi + q2)T^,^^\ql;m,m), etc. In the above 
equation, the magnetization m is taken as uniform 
(and the system is then uniform and isotropic at large 
distance) . The constant A is the variance of the random 
field, possibly renormalized to include fluctuations from 
the microscopic scale down to the beginning of the 
asymptotic regime. Identities similar to Eq. (j66|) are 
also derived for the higher orders. 

For the fiuid in the presence of a disordered porous ma- 
terial, Ward-Takahashi identities analogous to those for 
the RFIM should apply in the asymptotic regime, pro- 
vided of course the supersymmetry is not spontaneously 
brokenSi. For instance, when fc — >■ and g — >■ (and if 
fi < 2r]), one expects that, similarly to Eq. (|55| . super- 
symmetry implies that the disconnected direct pair cor- 
relation function (whose definition in the formalism with 
explicit replica-symmetry breaking will be made precise 
below) satisfies 



Cdrs,k{q^;p,p;po) 



K dq2 Cc 



,fc(9^;p;po), 



(67) 



with K a constant > which is essentially the quantity 
Cdis,k{q) - PoSmmiq)cmf,kiq)'^ evaluated for g = and 
renormalized only down to some k at the beginning of 
the asymptotic regime. 

Finally, in order to beyond the OZ/LPA approximation 
and introduce a renormalization of the field that allows 
us to study criticality in three dimensions, it is more 
convenient to replace the hard IR cutoff by a smooth 
one^S. Approximation schemes then exist for improving 
the description of the long-distance spatial dependence 
of the correlation functions and provide accurate values 
of the anomalous dimensions ?/ and fj—'—. 



VI. TOWARD A PROPER HRT FOR FLUIDS IN 
DISORDERED MEDIA 



A. Explicit replica-symmetry breaking 

To adapt the formalism with an explicit breaking of 
the replica symmetry (see above)^, we start with the 
grand partition function(al) of the fluid in the presence 
of a given configuration of the solid matrix characterized 
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by the microscopic density pq: 
Elfi-po] = exp{W[p,;pQ]) 
^J2^ [ - [ exp{-/3yfl({x,}^;po)+ 



N^ / 



(68) 



'/'(|2;"y|)p(a;)p(y)+/3 / p,{x)p{x)) 



X Jy 



where as before VR{{xai}N]Po) contains the sum of aU 
pairwise steep repulsive fluid-fluid interaction vr among 
atoms of replica a as well as the sum of all pairwise in- 
teractions Vmf between the fluid replica a and the solid 
matrix, and = —jSwff. 

We next consider n copies or replicas of the fluid in the 
presence of the same matrix configuration but coupled to 
different chemical potentials and we define 



expM^[{Ma};Po] = ]^^[Aia;po] == exp I ^W[p.a]Po] j 

a=l \a=l I 

(69) 
where the overline denotes an average over the disorder, 
i.e. the configurations of matrix particles po; po = pq is 
the average solid density and pa is the chemical poten- 
tial for replica a. The functional M^[{/ia}; po] can be ex- 
panded in cumulants of the random functional M^[/i;po]: 

M^[{Ma};H- 

n -. n 



,6=1 



(70) 



^ X! ^[Ma;po]M^[/ib;/5o]M^[/ic;po]| 



a,6,c=l 



A convenient trick to extract the cumulants of 
W[/Za;/5o] with their full functional dependence (i.e., 
generically, with different arguments pa, M6, etc.) is then 
to let the number of replicas be arbitrary and to view 
the expansion in the right-hand side of Eq. (170)) as an ex- 
pansion in increasing number of unconstrained, or "free" , 
sums over replica o^^i^^'^^ : 



n _, n 



Aifc;poJ 



a.h=\ 



1 " 
— ^ W■i{pa,Ph■.[i.c\PQ\^ 



3! 

a,&,c— 1 

(71) 
where the term of order p in the expansion is a sum over 
p replica indices of a functional Wp depending exactly 
on p replica sources, this functional being precisely equal 
here to the pth cumulant of W\p\pQ\. This procedure, 
which rests on an explicit breaking of the replica sym- 
metry through the introduction of the chemical poten- 
tials p,a is a 'priori different from the standard use of 
replicas, in which all chemical potentials (or sources in 
field-theoretical language) are equal; it therefore avoids 



the delicate handling of a spontaneous replica symmetry 
breaking when present'''^. The first cumulant W\ gives 
the grand potential of the fluid, W\ \p\ po] ~ — /3n/[/x; po] 
and the higher-order cumulants provide information on 
the probability distribution of the random functional 

w\p-rp^\ 

An effective action or Helmholtz free-energy functional 
is as usual defined through a Legendre transform. 



r[{pQ};po] = -w^[{Aia};po] +^X! / Po.ix)Pa{x) 

„ 1 J X 



where 



Pa{x) 



m{{paYp^\ 

PSpaix) 



(72) 
(73) 



is the fluid-density fleld in replica a. (As before, if neces- 
sary, one can take into account self-interaction terms by 
replacing pa{x) by pa{x) - Wff{Q)/2.) 

Just as the functional W, the effective action can be 
expanded in increasing number of free replica sums, 

n -. n 

r[{pa};po] =ro[po] + ^ri[pa;po] - - Yl r2[pa,pb;po] 



,.fc=l 



1 ^ r 

^ Z^ ^3[Pa,Pb,Pc 



3! 



Pol 



a,6,c— 1 



(74) 



where for later convenience we have introduced a minus 
sign for all even terms of the expansion. We have added 
a first term, Fq, that only depends on the solid matrix 
and is useful to generate the direct correlation functions 
(IPI vertices) of the matri»2L. 

r[{pa} ; po] and Ty[{pa} ; Po] are related by a Legendre 
transform, so if one also expands the sources Ha[{pf} ; po] 
(where we have denoted {p/} the n replica fields to avoid 
confusion in the indices) in increasing number of free 
replica sums, one can relate the terms of the expan- 
sion of the effective action to the cumulants of the ran- 
dom functional W^[J;/i]. The relation is straightforward 
for the first terms, but gets more involved as the order 
mcreases=^. 

More precisely, ri[p;po] is the Legendre transform of 
Wi[J;po], namely. 



with 



ri[p;Po] = -Wi[p;pq] + / p{x)p{x), 



p{x) 



SWi[p;pq] 



Sp{x) 
and the second-order terms is given by 

r2[pi,P2;po] = W2[p[pi; Po], p[p2; Po]; Po], 



(75) 



(76) 



(77) 
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where fj,[p; po] is the nonrandom chemical potential de- 
fined via the inverse of the Legendre transform relation 
in Eq. ^^, i.e., 



Legendre transform 



^i[p]Po]{x) = 



5p{x) 



(78) 



(Note that p[pa', po]{x) introduced here differs from the 
som-ce /l<q(x) appearing in Eq. (17^ : through the Legendre 
relations, the latter depends on all the replica fields {p/} 
while the former depends on a single replica field Pa-) 
The above expression motivates our choice of signs for the 
terms of the expansion in free replica sums of r[{pQ} ; po]i 
Eq. (1741) : r2 [pi , p2 ; po] is directly the second cumulant of 
W[ii] po] (with the proper choice of p[p; po]). 

We point out that rp[pi, ..., pp] for p > 3 cannot be di- 
rectly taken as the pth cumulant of a physically accessi- 
ble random functional, in particular not of the disorder- 
dependent Legendre transform of W[p\ h] (although it 
can certainly be expressed in terms of such cumulants of 
order equal to, or lower than, p). In the following and by 
an abuse of language, we will nonetheless casually refer 
to the Fp's as "cumulants of the renormalized disorder" 
(which is true for p = 2). 



B. HRT w^ith an explicit breaking of replica 
symmetry 

We are now in a position to formulate a more power- 
ful HRT framework for the quenched-annealed mixture. 
This first goes through the definition of a grand partition 
function(al) at the scale k, 

Sfc[{Ma};Po] =exp(VFfe[{pa};/5o]) = Yl AT \ M \ 



[ ... / exp{-/3^1/H(KJ 



N^;po)+ 



n 



a.b=i-"'-'y 



\X - V\)5ab - (l>k,ab{\x - V\)\pa{x)pb{y) 






(79) 

where 4>k,abil) = ~Pwk.,ab with Wk,ab{q) == 
'Ti-k.abiq)wff{q), a,b = l,...,n. Note that we do 
not consider any cutoff function for the matrix-matrix 
and matrix-fluid interactions. 

We next introduce the generating functional of the 
Green's functions at scale fc, W^fcij/ia}; Po]i through 



expW^fe[{^a};Po] ^ ^k[{Pa}; Po] ^C^p{Wk[{pa};Po]) 

(80) 
and the effective average action (Helmholtz free-energy 
functional) at scale k, Tk[{pa}', Po], through the modified 



n „ 
^k[{Pa};Po] = ~Wk[{pa}; Po] + PYI / Paix)paix) 

+ 9$! / M'k,abi\^~y\)Paix)Pb{y)- 



a.b=l''^-'y 

(81) 

The evolution with k of the effective average action 
again follows an exact fiow equation. 



dt^k[{Pa};Po] = 



9 X! / I ^t'i^k,ab{\x-y\)Fk,ab[x,y;{pa}]Po\, 



(82) 



a,fc=i-'^-'y 



where the 2-point (Green's) correlation matrix at scale 
k, Fk = W^ , is the inverse in the sense of matrices and 

operators of F^Jx,?/; {pa};Po] + (t)k;ab{\x - y\). 

The reasoning developed in the previous subsection can 
be applied to the effective average action Tk and its ex- 
pansion in free replica sums. As a results, Eqs. ([70]) to 
(|78p can be extended to any running scale k. To make 
the expansion in free replica sums a fully operational pro- 
cedure, one needs to be able to perform systematic alge- 
braic manipulations, as for instance the inversion of the 
matrix F^ +4'k- We detail in appendixiBJthe method for 
matrices that depend on two replica indices but are func- 
tionals of the n replica fields; extension to higher-order 
tensors is presented in Ref. |69j . 

From the detour via the field theory of the RFIM and 
its supersymmetric formulatioi>2i, we know that, con- 
trary to what was done in the naive HRT approach, 
the matrix (jjk^ab cannot be simply taken as diagonal 
despite the fact that there is no direct interaction be- 
tween atoms belonging to different replicas. The repli- 
cas are indeed indirectly correlated through the inter- 
action with the quenched matrix and the corresponding 
fiuctuations must also be regularized. We have seen in 
section IV Bl that this is necessary to avoid an explicit 
breaking of the underlying supersymmetry. We therefore 
choose (l)k,ab{q) = -l3Rk,ab{q)wff{q) with 



Rk,ab{q) = Rk{q)5ab + Rk{q) 



(83) 



with Rk{q) oc —dq2Rk[q). We shall denote by 0^ and 0fc 
the corresponding elements of the matrix (pi^ . 

Using the results of the previous section and of ap- 
pendix [B1 we derive through systematic expansions in 
free sums over replicas a hierarchy of exact HRT flow 
equations for the "cumulants" associated with the effec- 
tive average action (Helmholtz free-energy functional). 
The functional equation for the flrst cumulant reads 



dt^ki [pwpo] = 2 



dt{<i>k{q) + 4'k{q))Fcon.k [pi;pa] 



-qq 



dt<j>k{q)Fdis,k [pi,pi;po]_ 



(84) 
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where the correlation functions (or Green's functions or 

else "propagators") Fcon,k = Fk and Fdis,k = Fk are 
obtained as zeroth-order components of the expansion in 
number of free replica sums of the matrix Fk] they are 
explicitly given by 



the quenched-annealed mixture in the explicit replica- 
symmetry breaking framework. We consider for simplic- 
ity uniform replica density fields, which leads to 



/ (2) ^ 

FconAPii Po] = {J^ki [Pi'^Po] + 4>y 



FdisAPl^P'i'iPQ] = ^con,fc[Pl] I r[,2 [Pl, /O2; Po] 



(85) 



/OqX 



.(i;i) 



r^l [pi] Po\Sm,n[po\T\{ '[P2]P0\ 



(l;l)r 



Fcon,k[P2': Po], 

(86) 

where all quantities have to be considered as operators 
in real or Fourier space when the density fields are inho- 
mogeneous. The above expressions follow directly from 
the results derived in appendix [BJ Note that the 2-point 
IPI vertex fJ^j^' [pi] po] only involves the first, 1-replica, 
cumulant and corresponds to a matrix-fiuid direct corre- 
lation function in the language used in sections IIIII and 

lYl 

The functional HRT equation describing the flow of 
the second cumulant can be expressed as 



<9tFfc2 [pi,p2;pa] = T;9tTr\ Fk [pi] Po] (r^j"' [pi, P2; Po] 



(110) 



-.(20) 



fc3 [pi,Pi,P2;po]) +Fk[pi,pi;po]ri^' [pi,P2;po] 



.(11) 



(i;i)r 



■^Fk [pi,p2;po] (r^,2 [pi.P2;po] + po^i{ '[pi;po]Smm[po]x 
rii^^[p2;po] - (fk) + y^fe [pi,P2;po]r^i^^[pi;po]S'mm[po] 



^(10;1) 
fc2 



(FL •[pi,P2;poj 



POj.(oa) 



(i;i)r 



perm{l2) 



^kl [P2; Po]Smni[Po]'^kl iP^', Po]) 



(87) 



where perm{12) denotes the expression obtained by per- 
muting pi and p2 (some care is needed in the term 
by term identification in order to properly symmetrize 
the expressions and satisfy the permutational property 
of the various arguments of the cumulants), the op- 
erator dt acts only on the cutoff functions (with e.g., 
dtFk-qq' = - J^„ Fk-q"dtRkiq")Fk-q"q'), aud the trace Tr 
is over momenta. Similar equations can be derived for 
the higher-order cumulants. 

As in the case of the naive HRT, the hierarchy for the 
cumulants of the renormalized disorder, when expressed 
in a functional form as above, contains full information 
on the complete set of IPI (or direct) correlation func- 
tions. Flow equations for the latter are simply obtained 
by taking appropriate functional derivatives. 

We conclude this subsection by providing the coun- 
terpart of Eqs ([551) ^^^ ((55)) . i.e., the OZ equations for 



1 + pi/icon,fc(g;pi;po) 



1 



1 - piCcon,k{q] pi; Po) 



hdis,k{q\Pl,P2]P(i) = [l + PlhconMi; Pi; Po)] X 

[1 + P2^con,fe(g;p2;po)] [cd2s,kiq;pi,P2;po) + PoSmmiq; po) 

X Cmf,kiq; Pi; pa)cr,if.kiq; P2; Po)], 

(89) 
whereas, the matrix fluid total pair correlation function 
is given by 



hrnf,k{q;pi;po) = Sr,im{q;po) 



Cmf,k{q; Pi; Po) 
1 - PiCcon,k{q; pi; Po) 



(90) 
To make contact with the replica-symmetric formu- 
lation [see Eq. (PB| ]. we have defined in the above 

(2) ^ 

expressions Ccon,k{l'^ P'^ Po) = ^kiil'P'Po) + (t>k{q)-, 

CmfAT,p;po) = '^ki^\q;p;po), and cdis,k{^; pi, P2; po) = 

FJ.2 (q; Pi,P2; Po) — 't'kiq)- The disconnected direct cor- 
relation function Cdis.k{q; p) of the replica-symmetric OZ 
equations is equal to Cdis,k{q;pi = P,P2 = p; Po) in the 
above formalism; the correspondance is obvious for the 
two other direct correlation functions. 



C. Asymptotic analysis 

As in the case of the naive HRT, the hierarchy simpli- 
fies when considered in the asymptotic regime fc,q — > 
near the critical point. We choose the IR regulators such 
that 



Rk{q) = Zkk'Hr) 



and 



^f^.2\ 



Rk{q) = Kkf{q 



-(2')-')) 



(91) 



(92) 



with Zk - k^-^, Kk - fc~ 

Similarly to what was done before, we introduce di- 
mcnsionless quantities (after accounting for scaling di- 
mensions that are appropriate near a zero-temperature 
fixed point): ipa ^ k-'^'^-'^+'^^^ipa - Pcrit), afe(<p) ~ 
/c-('^-«)[yl/,fe(p; Po) - Af^k;crtt{p; Po)] with e = 2 + ij-fi, 
etc. [We recall that Af.k{p;PQ) = Tki{p; po)/V.] We as- 
sume for simplicity in the following that 77 < 2r7 so that 
all contributions involving matrix-fluid correlations are 
subdominant and can be neglected when k ^>- 0. 

For a uniform density field, the first equation of the 
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HRT hierarchy becomes 



1 



dtakiif) + {d-2 - 7] + ri)akif) - -(d - 4 + ri)ifd^ak{(p) 



[Ucon,k{Sp]V) +f{q^)f 



1 



[ucon,fe(g^;<^) + ^('7^)] 



(93) 

where, to make a direct contact with the previous rephca- 
symmetric treatment (see section HV)) . we have defined 



„,(g2;^)^fc-(2-")rfi)(g;p;po) (94) 



and 



UdzsAf'^'Pi^f2)^k'^'^ T[."'((7;pi,/92;Po), (95) 

with the dimensionless fields tpa and the densities pa re- 
lated as above and q = kq. (Note that we have chosen 
to not include the regulator f in Udis,k to keep the de- 
pendence on f explicit, despite the fact that as defined 
below Eq. (^0)) . it is present in Cdis.k-) For compactness, 
we have also introduced the short-hand notations 



x'2''/^2\ 



and 



rit) = i2-v)fir)-2rr'{r) 



fif) = -i2r^-fj)rif)-2fP'{f). 



(96) 



(97) 



Note that all constants (that do not depend on k) can 
be incorporated in a trivial redefinition of the fields and 
functions. 

Along the same lines, one obtains the flow of 

vk{'fu^2) - fc-(''-2»)[rfc2(pi,P2;po)/v^] as 

dtVk{ipi,(p2) + {d-4-2T] + 2fj)vk{(pi,(p2) - -^{d - 4 + fy) 

X (vJi^^i +'p2d,^^)vk{(pi,ip2) 
1 f r r(g2) 



whereas Udis,k{q = 0; v^i, v'2) = d^^d^k{^i,V2)- 

From the discussion in section IV Bl one knows that 
the underlying asymptotic supersymmetry implies Ward- 
Takahashi identities. In particular, for the the IR cutoff 
functions the identity in Eq. (157]) implies that R{q) = 
—Kdq2R{q) with K a constant. We therefore take the 
dimensionless cutoff functions such that 



(110)/ -2 N 

-7fc3 (9 ;</'i>¥'i,¥'2) 



2 7g I [UconAi^'^ <Pl) + ^(9^)] 
,„ (20). .2 J"dis,fe(g^;<y5l,'^l) -r(g2)] 

[Mrfis,fc(g^;'jgi,¥^2) -r{q^)f 

[UconA'PWl) + 'r{¥')][ucon,k{q^]^2) + f{q^)] 

Hf) ( Ik2\i^'^fl^f2) 



(98) 



[Ucon,k{q^; ^l) + f{q^)] V [Ucon,k{q^]Vl) + f{Sp)] 

[Ucon,k[q]ip2) + r{q'')\ J 
where we have defined 

7if (<Z^;^i,^2) - fc^-'-'^rif (g;pi,p2;Po) (99) 

and 

7fc"°^(9^;¥'i>9'2,'P3) ~ fc^+^''"^''r{,""^(q;pi,p2,P3;Po), 

(100) 



Hr) = -dMr) = -r'(q^). 



(101) 



The supersymmetry is then satisfied when Kk = KZk- 
The above choice nonetheless allows for a spontaneous 
breaking of the supersymmetry when Kk ^ KZk- 

One can easily check that the above asymptotic HRT 
equations coincide with the nonperturbative RG flow 
equations for the RFIM22ri2i obtained in the explicit 
replica-symmetry breaking framework. As shown in the 
latter studies, the flow equation involving the dimension- 
less second cumulant, Eq. (|98p . turns out to be a key 
element to resolve the conundrum of dimensional reduc- 
tion. 



D. Supersymmetry-compatible OZ/LPA 
approximation 

We first reconsider the naive OZ/LPA approximation 
that we have used in the context of the replica-symmetric 
formalism (see section IIV Ep . Again, as there is no 
field renormalization and no disorder renormalization, 
r] — fj = 0, which implies that one cannot simply drop the 
matrix-fluid direct correlations in the asymptotic equa- 
tions. However, we assume that no accidental cancella- 
tion takes place and that Cdis,k{q)+ PoSmm{q)cmf,k{q)'^ — 
Kk — > Kq in the small-fc, small-g limit. We then repeat 
the derivation of section IIV El but we now consider two 
cutoff functions, Rk and Rk, which we can choose to be 
related so that their form does not explicitly break the 
underlying asymptotic supersymmetry. Accordingly, we 
use the Ward-Takahashi identity to require that Rk (q) = 
—Kdq2Rk{q), with K a constant to be determined later 
on. With the choice of regulator Rk{q) — k^f{q'^/k^), 
one therefore has Rk{q) = —Kf'{q^/k^). 

The asymptotic form of the HRT equation for the 
Helmholtz free-energy density ak{^) now reads: 

dtttkiv) + {d- 2)ak{(p) - x(d - 4)(pd^ak{(p) = 



Ko f 2[r{f)-ff'if)] K 

2 Jq[uconAf;^)+rm? 2"" 

2ff"{f) ^ 2[fie)-fr'if)]f'if) 

[Ucon,fc(g^; v) +?^(g2)]2 

(102) 



4 {uconAq'^'^'P) +riq^) 



with Ucon.kiq^) = q^ + a'kif). 

Let us introduce the variable y ^ q" and define as 
before Vd = [2''+^7r'^^^T{d/2)]-'^ [it is equal to l/4th 
of the area of the d-dimensional sphere of radius unity 
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divided by (27r)'']. Then, after using the fact that 
dy[ucon.k{y'' f) + ^(y)] = 1 + r'{y), one finds 



dtak{ip) + {d- 2)ak{(p) - -{d - 4)ipd^akiip) 



VdiKo- K) dyy2 



VdK / dyy-' ^dy 



A-i 2[f(y) - 2/f'(y)] 



2[f{y)-yf'{y)] 



[uconAy^v) + f{y)] 



After integration by parts, the second term of the right- 
hand side can be reexpressed as 



VdK dyyi ^d, 
Kvd 



2[f{y) - yr'(y)] 
^ '' 'uco„,k{y;^) + f{y)]^ 
^-1 '2[r{y)~yr'{y)] 



dyy~ 



[ucon.kiy'^'p) + fiy)] 



\4nj2j^ Uco«,fc(g2;(^)+f((72)' 

(104) 
where the integral over the wavevector q is now in a (d — 
2)-dimensional space. In deriving the above equahties we 
have used the fact that Vd{d — 2)/2 = Wd-2/(47r). 

If we now choose K — Kq, the HRT equation for ak{(p) 
reduces to 



dtak{ip) + (d - 2)ak{(p) - -{d 



Ka 

in 



1 f'^''-^^ 2[r{f) 






(105) 



e 



<{^) 



f{q^y 



which, up to the trivial constant factor Ko/{An), is 
exactly the HRT equation for a pure fluid in the 
OZ/LPA approximation in dimension d—2 [compare with 
Eq. (1581) ]. Within this approximation, the exponents de- 
scribing the gas-liquid critical point of a fluid in a disor- 
dered porous material are therefore identical to those of 
the pure fluid in two dimensions less. For instance, for 
the sharp cutoff in d — 5, this provides the exponents 
V ~ 0.689, 7 ~ 1.378, 5 = 5, and tj ~ 0.581 as obtained 
for the bulk fluid in d = 3 in the same approximation and 
with the same sharp cutoff^ '2. This proves that dimen- 
sional reduction is obeyed within this OZ/LPA approxi- 
mation, provided that one properly chooses the infrared 
regulators. The results from the naive OZ/LPA approx- 
imation were thus plainly wrong. 



E. Improved truncations 

The OZ/LPA approximation, even when designed not 
to break the supersymmetry explicitly, is not satisfac- 
tory. First, it does not allow one to investigate a spon- 
taneous breaking of supersymmetry and the associated 
breakdown of dimensional reduction at the critical point, 
and second, it prevents one from studying dimensions 



less than 4. Improved truncations of the HRT hierar- 
chy are therefore required which at least include the sec- 
ond cumulant of the renormalized disorder for arguments 
that are generically different (so that a nonanalyticity can 
freely emerge along the RG flow, see above) and provide 
a better description of the long-distance behavior of the 
correlations. 

From the nonperturbative RG work on the RFIM22ri2i, 
we have learned that a truncation must include the 
first cumulant r^i [p; po] and the second cumulant 
^k2[piT P2'i Po\- In the HRT framework, which is oper- 
ationally based on coupled fiow equations for IPI vertex 
functions evaluated for uniform density fields, this im- 
plies approximation schemes to the two-point (or higher- 
order) direct correlation functions obtained from T^i and 
Tk2- Interestingly, the formalism with explicit breaking 
of the replica symmetry that we now use provides us 
with additional "compressibility sum rules" . At the level 
of the direct pair correlation functions, in addition to 
Eqs. ([57)1 and ([55]) relating the connected fluid-fluid pair 
correlations and the matrix-fluid ones to the flrst cumu- 
lant, i.e., the Helmoltz free-energy density of the fluid 
AfAP'iPo) = ^ki{p\ Po)lV , there is a new sum rule for 
the disconnected fluid-fluid direct pair correlation func- 
tion. 



Cdis,k{q = 0; pi, P2; Po) ~4>k{q = 0) 



= r 



(11) 

fc2 



{pi,P2;po) 



d'^Vkipi, P2] Po) 



(106) 



dpidp2 

where Vk{pi, P2; Po) = '['k2{pi, P2; Po)/V is the second cu- 
mulant for uniform fields (divided by the sample volume 
V whose notation should not be confused with Vk). 

A minimal truncation of the HRT hierarchy should 
thus include Af^k{pi;Po), Vk{pi, P2; Po), Ccon,k{'i', Pi' Po), 
Cmf,kil',Pi',Po), and Cdis,k{<l', Pi, P^', Po), with sum rules 
relating the three direct correlation functions at g = 
to the two potentials. The remaining challenge is to 
go beyond the OZ/LPA treatment of the spatial or q 
dependence of the correlation functions. In the field- 
theoretical framework, there are two ways to achieve this 
and, in consequence, to predict nonzero anomalous di- 
mension(s). The standard procedure is to consider next 
orders of the "derivative expansion" , which is an expan- 
sion of the effective average action in gradients of the 
field^. A recently implemented alternative consists of 
truncating the hierarchy of equation for the IPI ver- 
tices by approximating the momentum dependence of the 
3- and 4-point vertices so that the equation for the 2- 
point vertex can be closed through a systematic use of 
compressiblity-like sum rules^^. (This is akin to a proce- 
dure formulated earlier by Parola and Reatto^ for fluids 
but never implemented in this context.) In the context 
of the RFIM, it has however proven difficult to use this 
last approximation as it explicitly breaks the underlying 
supersymmetry^s. 

It appears that a possible candidate for a practical 
implementation of the HRT for describing fiuids in dis- 
ordered porous materials from short- to long-distance 
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physics would be a combination of the derivative expan- 
sion for small wavevectors with an OZ-like approxima- 
tion carrying all microscopic details for high wavevec- 
tors. Dropping for simplicity the explicit po dependence 
in the direct correlation functions, a typical approxima- 
tion would then read 



1 
P 

d^Af,k{p;Po) 



At,p) ~ (t>k{q) 



5p2 



+ Zfc(p; po)q^Q{q - q^) - [1 - e(g - q^)] X 



[ccon,R{q) - Ccon.Riq ^ 0) + 4>{q) - (t>{q = o)] , 

(107) 
where Q{q — q^) is a smoothed Heaviside step function 
that goes continuously from 1 when q < q^ to Q when 
q > q^ with g^ a crossover wavevector larger than the IR 
cutoff k and smaller than the typical interatomic distance 
in the liquid, a [one could also use an approximate form 
in the spirit of Eq. (fTS)) ]. and 



Crnf,k{q;p) = - 



d'^Af^kip;Pa) 
dpdpo 



(108) 



+ [CmfMq'^ P) - Cmf,R{q = 0; p)], 



Cdis,k{.q\ Pi, P2) - 4>k{q) 



d^Vk{puP2) 



dpidp2 

Cdis,R{q = 0;pi,p2) 



+ [cdis.R{q\PliP2) 



(109) 

Again, one can improve this ansatz by implementing a 
core condition for the total fluid-fluid and the matrix- 
fluid pair correlations. In addition, we set to zero the 
third and the higher-order cumulants, so that the HRT 
hierarchy is closed. 

In the asymptotic regime, when fc, g ^- (and if 
77 < 277 so that one can drop the matrix-fluid contribu- 
tion), one has Ucon,fc (9) ^ a'l.{ip)+Zk{^)q^ a.T\dudis,k{q) = 
dipidip^Vk{(pi,^2)- It is easy to show^i that the super- 
symmetry is satisfied when d^^d^^Vk['^iT^2)\^i=^2=v = 
Zk{y^) and (spontaneously) broken otherwise. This trun- 
cation has been shown to provide a very good description 
of the critical behavior of the RFIM in d = 3 while re- 
solving the dimensional- reduction puzzle^'*. 

With the above closure, one has to solve, on top of 
the OZ equations, three coupled nonlinear partial dif- 
ferential equations describing the flows of Af^kip', Po), 
Vk{pi, P2', Pq) and Zk{p',po), the latter being obtained 
from dq2T f.-^ (q-^; p; pQ)\q^Q. The solution would allow a 
complete description, at all length scales, of a fluid in 
a disordered porous medium in rf = 3. However, this 
represents a very arduous numerical task2^. 



VII. CONCLUSION 

In this article we have considered the equilibrium be- 
havior of fluids adsorbed in disordered mesoporous ma- 
terials, with a special attention to the gas-liquid critical 



point. Because it allows one to keep track of the details 
of the physics at a microscopic level while dealing with 
fluctuations on arbitrarily large length scales, the HRT 
provides an efficient framework. A nontrivial general- 
ization of the HRT of bulk fluids and mixtures is how- 
ever required. Through a description of the fluid/matrix 
system as a quenched-annealed mixture, we have com- 
bined liquid-state statistical mechanics and formalism 
borrowed from the theory of systems with quenched dis- 
order. 

A straightforward implementation of the HRT to the 
replica description of the quenched-annealed mixture, 
which in particular assumes replica symmetry, has been 
shown to lead to results concerning the critical behav- 
ior that are inconsistent and unsatisfactory. At the same 
time, it has provided strong indication that the critical 
behavior of the quenched-annealed mixture and that of 
the RFIM are in the same universality class. The same 
difficulties encountered in the RG treatment of the lat- 
ter are then to be expected in the HRT of the former. 
We have then built on the recent nonperturbative RG 
approach of the RFIM developed by two of us for solv- 
ing pending puzzles in the long-distance behavior of the 
mode]32ri2i to propose a more sophisticated HRT for the 
quenched-annealed mixture. As discussed and illustrated 
in the paper, this formalism opens the way to cure the 
inconsistencies of the naive treatment. We have also de- 
vised approximations to the HRT hierarchy that are ex- 
pected to lead to an accurate description of the behavior 
of fluids in disordered porous materials. A full-blown 
resolution of the resulting set of coupled nonlinear dif- 
ferential equations is however a formidable task that we 
defer to future work. 



Appendix A: Higher orders of the naive HRT 
hierarchy in the asymptotic regime 



We have derived in section IIVCI the asymptotic equa- 
tion for the Helmholtz free-energy density and we have 
introduced scaling dimensions appropriate for studying a 
zero-temperature fixed point. The scaling of the 2-point 
correlation function functions is given in Eqs. (I40p and 
(|45)) . Here, as often in the main text, we assume for 
simplicity that fj < 2?], which implies that the contri- 
bution from the matrix-fluid correlations always lead to 
subdominant terms in the asymptotic regime. 

Some care is needed to derive the proper scaling of the 
higher-order direct correlation functions as there are a 
priori an increasing number of distinct functions when 
n -^ (whereas only the connected and disconnected 
components appear at the pair level). For instance, the 
3-body direct correlation functions should scale as 



i™o IZrfeS(9i'«2.93) 



yijy2,y3; 
.,b=i (Al) 

k''-'-i^''-^+^^u,Aq,,q2,h) 
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I^1^0^^lll2a(^^'^^''i3) 



(A2) 



Similarly, 

n 



(A3) 



n^O 



a.b,c=l (A4) 



i^n XI ^M2ab (91 , 92 , 93 , 94) 



"^"a^l '^'"^ (A5) 

^-^'^-^'-^('^-4+'^) ^12,^(91,92, 93, 94), 

etc., with 9 = 2 + T] ~ fj. The expressions of the various 
functions appearing in the above equations can actuahy 
be rationalized by making use of the cumulants of the 
renormalized disorder that can be introduced when ex- 
pUcitly breaking the rephca symmetry (see section IVI Al 
and the following appendix). This connection also justi- 
fies the choice of sign used in the above equations. 

We illustrate the higher-order equations of the naive 
HRT hierarchy for case of the connected direct pair cor- 
relation function (2-point IPI vertex). It reads 



dtUcon,k {q^\ip) + {2- 1])UconM (9^ ; <y3) - - (d - 4 + ?]) X 



d 



2 V^a [uconAq'^)+r{q'^)? 

2 K,fc(9, q', -~{q + q'))? u^sAq'^) 

[uconMq'^) + riq'^)?[uconM\q + 9'P) + ri\q + 9?)] 

g Mi2,fc(9', -(9 + q'),q) ui.kjq', -(9 + q'),q) 

[uconAq'^) + r{q''^)][uconA\q + 9T) + ri\q + 9T)] 
Wi2,fc(9',-9',9, -9) 



+ 



[uconAq'^) + r{q")]j' 



(A6) 

where the dependence on the dimensionless field ip is not 
explicitly shown in the right-hand side and we recall that 
the operator dt only acts on the regulator, with the for- 
mal definition dtr{q^) = r{q^) = (2 — ri)r{q^) — 2q^r'{q^). 
Similar equations can also be derived for the fiow of 
Udis.k{q^','p) and all higher-order dimensionless IPI ver- 
tices but it is not worth displaying them here. 

The existence of "compressibility sum rules" and their 
generalizations, which merely express the fact that the 



effective action is the generating functional of the IPI 
vertices, allows one to identify the above equation for 
Ucon.kiq'^]^) when 9^ = with the second derivative 
with respect to the field ip of the flow equation for ak{(p), 
Eq. (gT]). With the help of the operator dt, Eq. (gT]) can 
be reexpressed as 

dtUkiip) + {d-2 - ij + ri)ak{(p) - -{d - A + fj)ipd^akiip) 
Udis.kiq"^) 



[ucon.kiq^) +^(9^)]' 

' (A7) 

Then, by deriving twice this equation with respect to ip, 
one obtains 

dta'^iip) + {2- r,)al{^) -\{d-A + fj)ipd^4{v) = 



¥' 



udis,k{q^)cil,ucon,k{q^) d'iudisAq^) 



.,km + r{q^)] 



[Ucon.k{q^) +^(9^)] 



UdtsAq^)[^v'^con.k{q^)? 



, d^udisAq^)(^v>'^con.k{q'^) 

[Ucon,fc(9^) +?'(92)]2 

where again the if dependence is omitted in the 
right-hand side. The compressibility sum rule for 
the connected fluid-fluid pair correlations states that 
UconAf = 0;<^) = d^akltf) = a'^if)- Then, by set- 
ting 9^ = in Eq. (|A6[) and using the generalized sum 
rules, ui,fe(0,9, -9) = d^UconA'f')^ wi2,fe(9, -9, 0) == 
d,^Udis.k{q^)^ wi,fc(9,-9, 0,0) = 9^Wcon,fe(9^), and 
wi2,k{q^ -9' 0, 0) = dludisAq^)^ o^e recovers Eq. (JA8|- 



Appendix B: Algebraic manipulations for the 
expansions in number of free replica sums 

We consider the replicated version of the quenched- 
annealed mixture in which the chemical potentials are 
different for each fluid replica so that the density flelds 
{pa] are also different. The effective action (Helmholtz 
free-energy functional) for the (n -I- l)-component mix- 
ture can be expanded in increasing number of free replica 
sums as indicated in Eq. (TMl) . The developments below 
generalize to the (n-l-l)-component system the results for 
the RFIM—; in the latter the average over the random 
field has been performed so that one is left with an n- 
component mixture with, however, interactions between 
the replicas. 

Consider a generic (symmetric) matrix Aap[{p f} ] po], 
where we have denoted {p/} the n replica density fields 
to avoid confusion in the indices (we recall that Greek 
indices denote any component, including the matrix, 
whereas Roman indices denote fluid replicas only). Its 
nx n submatrix associated with the fluid replicas can be 
decomposed as 

Aab[{pf};po] ^ Aa[{pf} ; po]Sab + Aab[{pf} ; po]- (Bl) 

In the above expression, it is understood that the sec- 
ond term Aab no longer contains any Kronecker symbol. 



21 



Each component of the matrix A can now be expanded 
in increasing number of free rephca sums, 

n 

^oo[{p/} ; Po] = 4"(] [Po] + E 4'cl be; Po] + • • • (B2) 



^Oa[{P/};Po] = A[[j'[Pq;Po] +X!'^Ol[Pa|Pc;Po] 



c=l 



(B3) 



X[{p/};po] - A[°l[p,;po]+E^™[^''l^-;^o] + --- (B4) 



c=l 



c=l 

(B5) 
where the superscripts in square brackets denote the or- 
der in the expansion (and should not be confused with su- 
perscripts in parentheses indicating partial derivatives). 
As an illustration, the expansion of the matrix FJ. of 
2-point IPI vertices (direct correlation functions) reads 
in terms of the expansion of effective average action it- 
self (for ease of notation we drop the subscript fc in the 
expressions appearing in the right-hand sides): 



rSo[{p/};po] = r(^^[po] + Err^[/'-^o] + --- (B6) 



ri;Oa[{P/} ; Po] =- ^x [Pa; Po] "Xl ^2 ' \Po.. Pc\ Po] 



ri;l[{p/};po] = rf^[pa;po] -Er^^°^[p„,pc;po] 



(B7) 



(B8) 



^(2) 



(ll)r 



n 

+ H^3"°^[Pa>P6>Pc;Po] H , 

c=l 

(B9) 

where the permutational symmetry of the arguments of 
the Ffep's has been used and the superscripts indicate the 
number of functional derivatives with respect to the den- 
sity fields appearing as arguments. For simplicity, we 
have not used the superscript when there is no deriva- 
tive with respect to po and we have not indicated the 
dependence on spatial coordinates or momenta. 

The above expressions can be used to reinterpret the 
various direct correlation functions (IPI vertices) appear- 
ing in the replica-symmetric formalism. We focus here 



on the fluid-fluid correlations, extension to the matrix- 
fluid ones being straightforward. One readily obtains 
that Ccon,ki<l, p) is the same in the two formalisms and 
is related to a 2-point IPI vertex obtained from the first 

(2) 

cumulant, Tl.^{q;p); Cdis,k{Q',p) in the replica-symmetric 
formalism is equal to Cdis,k(<l', Pi = P, P2 = p), which is 
related to a 2-point IPI vertex obtained from the second 

cumulant, F^j ('ZSPiP); in the formalism with explicit 
replica symmetry breaking (for Ccon,k and Cdis,k one also 
has to include IR regulators). 

Similarly, the higher-order vertices appearing in the 
previous appendix can be expressed as 



lim 



^-^ 

a,b=l 



rl;fl6^(9i,92,93;p) = ^kii<ii:i-2ii3;p), (Bio) 



!™nXlrM2a^(^l'92,93;p) 



n^O 



(Bll) 



('21 ^ f 1 1) 

-[^k2 (9i.92,g3;p,p) + r^2 (9i>92,g3;p,p)], 



iim,rla^2l^(9i,92,g3;p) = r[.3"^(gi, 92, 93; p,p,p), 

(B12) 
etc. 

Algebraic manipulations on matrices such as those de- 
fined above can be performed by a term-by-term iden- 
tification of the orders of the expansions. For instance, 
the inverse B ~ A^^ of the matrix A can also be put 
in the form of Eq. (JB1[) and its components can be ex- 
panded in number of free replica sums. The term-by- 
term identification of the condition AB — 1 leads to a 
unique expression of the various orders, Bqq, Bq^, 

and B^P\ of the expansion of B in terms of the ^qq 

aI^I's, A[«1's and JlI'I's with q < p. The algebra becomes 
rapidly tedious, and we only illustrate here the results 
for the zeroth-order terms: 






Bol[Pl:Po] 



i3io'[Po]=4°J[Po] 



(B13) 



4o[po]<1pi;po]B[°i[pi;Po], (B14) 



m[p,;po]=A^%i;por\ (B15) 

sM [pi, P2; Po] = -5t"l [pi; Po] (^I[°l [pi, P2; po]- 

^fj [Pi; Po]4°J [Po]^S [P2; po]) b™ [P2; po]. 



(B16) 
Note that by combining Eq. (|B6p and Eq. (jBlSp . one 
can immediately see that for a uniform density field, 

^1^00 (9; Po) i^ simply l/[poS'mm(9;Po)] and is indepen- 
dent of the presence of the fluid. 
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